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Abstract. This paper is concerned with some nonUnear propagation phenomena for 
reaction-advection-diffusion equations with Kolmogrov-Petrovsky-Piskunov (KPP) type non- 
Hnearities in general periodic domains or in infinite cyhnders with oscillating boundaries. 
Having a variational formula for the minimal speed of propagation involving eigenvalue 
problems ( proved in Berestycki, Hamel and Nadirashvili [3]), we consider the minimal speed 
of propagation as a function of diffusion factors, reaction factors and periodicity parameters. 
There we study the limits, the asymptotic behaviors and the variations of the considered 
functions with respect to these parameters. Section 9 deals with homogenization problem as 
an application of the results in the previous sections in order to find the limit of the minimal 
speed when the periodicity cell is very small. 
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1 Introduction 

This paper is a continuation in the study of the propagation phenomena of pulsating travel- 
ling fronts in a periodic framework corresponding to reaction-advection-diffusion equations 
with heterogenous KPP (Kolmogrov, Petrovsky and Piskunov) nonlinearities. We will pre- 
cisely describe the heterogenous-periodic setting, recall the extended notion of pulsating 
travelling fronts, and then we move to announce the main results. Let us first recall some 
of the basic features of the homogenous KPP equations. 
Consider the Fisher-KPP equation: 

ut -Au = f{u) inM^. (1.1) 

It was introduced in the celebrated papers of Fisher (1937) and in [19] originally motivated 
by models in biology. Here, the main assumption is that / is, say, a function satisfying 

/(O) = /(1) = 0, /'(1)<0, f(0)>0, 

/>Oin(0,l), /<Oin(l,+oo), ^"'^^ 

/(s)</'(0)s,VsG [0,1]. (1.3) 

As examples of such nonlinearities, we have: /(s) = s(l — s) and /(s) = s(l — s^). 

The important feature in (1.1) is that this equation has a family of planar travelling 
fronts. These are solutions of the form 

V(t, x) G M X M^, u{t, x) = 0(x ■ e + ct), . . 

0(-oo) = and 0(+oo) = 1, ^ ' 

where e G M.^ is a fixed vector of unit norm which is the direction of propagation, and c > 
is the speed of the front. The function </) : M i— > M satisfies 

P{-oo) = and (f){+oo) = 1. ^ ^ 

In the original paper of Kolmogorov, Petrovsky and Piskunov, it was proved that, under the 
above assumptions, there is a threshold value c* = 2a//'(0) > for the speed c. Namely, no 
fronts exist for c < c*, and, for each c > c*, there is a unique front of the type (1.4-1.5). 
Uniqueness is up to shift in space or time variables. 
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Later, the homogenous setting was extended to a general heterogenous periodic one. 
The heterogenous character appeared both in the reaction-advection-diffusion equation and 
in the underlying domain. The general form of these equations is 

r Ut = V- {A{z)Vu) + q{z) ■ Vu + f{z, u), t z e 

\vA Vu{t, z) =0, t e R, z e on, ^ ^ 

where 1^(2;) is the unit outward normal on dQ at the point z. 

The propagation phenomena attached with equation (1.6) has been widely studied in 
many papers. Several properties of pulsating fronts in periodic media and their speed of 
propagation were given in several papers ( Berestycki, Hamel [2], Berestycki, Hamel, Nadi- 
rashvili [3], and Berestycki, Hamel, Roques [5, 6] and Xin [36]). In section 2, we will recall 
the periodic framework and some known results which motivate our study. The main results 
of this paper are presented in sections 3 to 6. 



2 The periodic framework 

2.1 Pulsating travelling fronts in periodic domains 

In this section, we introduce the general setting with the precise assumptions. Concerning 
the domain, let > 1 be the space dimension, and let d be an integer so that 1 < d < N. 
For an element z = (xi, X2, ■ ■ ■ , Xd, Xd+i, ■ ■ ■ , x^) G M^, we call x = (xi, X2, ■ ■ ■ , xa) and 
y = {xd+i, ■ ■ ■ ,X]\f) so that z = (x, y). Let Li, - ■ ■ .L^hed positive real numbers, and let VL 
be a nonempty connected open subset of satisfying 

{3i? > 0;V(x,?/) e 1], |?/| < i?, 
y{ki,---,kd)eLiZx---xLdZ, n = n + ^kid, ^'^'^^ 
k=l 

where {ei)i<i<N is the canonical basis of M^. In particular, since d > 1, the set Q is 
unbounded. 

In this periodic situation, we give the following definitions: 

Definition 2.1 (Periodicity cell) The set C = {{x,y) E Vl; xi G (0,Li),--- ,Xd G 

{0,Ld)} is called the periodicity cell ofQ. 

Definition 2.2 (L-periodic flows ) A field w : Q ^ is said to be L-periodic with 

respect to x if w{xi + ki, ■ ■ ■ ,Xd + kd,y) = w{xi, - ■ ■ ,Xd,y) almost everywhere in Q, and for 

d 

all k = {ki, ■ ■ ■ ,kd) E Y\ LiZ. 

1=1 

Before going further on, we point out that this framework includes several types of simpler 
geometrical configurations. The case of the whole space corresponds to d = N, where 
Li, . . . , Ltv are any positive numbers. The case of the whole space with a periodic array 
of holes can also be considered. The case d = 1 corresponds to domains which have only one 
unbounded dimension, namely infinite cylinders which may be straight or have oscillating 
periodic boundaries, and which may or may not have periodic holes. The case 2 < d < N — 1 
corresponds to infinite slabs. 
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We are concerned with propagation phenomena for the reaction-advection-diffusion equa- 
tion (1.6) set in the periodic domain Q. Such equations arise in combustion models for flame 
propagation (see [27], [31] and [37]), as well as in models in biology and for population dy- 
namics of a species (see [14], [18], [20] and [28]). These equations are used in modeling the 
propagation of a flame or of an epidemics in a periodic heterogenous medium. The passive 
quantity u typically stands for the temperature or a concentration which diffuses in a pe- 
riodic excitable medium. However, in some sections we will ignore the advection and deal 
only with reaction-diffusion equations. 

Let us now detail the assumptions concerning the coefficients in (1.6). First, the diffu- 
sion matrix A{x,y) = {Aij{x,y))i<ij<N is a symmetric C^'^{Q) (with 6 > 0) matrix field 
satisfying 



N 



A is L-periodic with respect to x, 
30 < ai < a2,V(x,y) G fi,V^ G 



(2.2) 



The boundary condition i/-AS7u{x, y) = stands for i^i{x, y)Aij{x, y)dx^u{t, x, y), and 

V stands for the unit outward normal on d^. We note that when A is the identity matrix, 
then this boundary condition reduces to the usual Neumann condition d^u = 0. 

The underlying advection q{x,y) = {qi{x,y), ■ ■ ■ ,qN{x,y)) is a C^'^{Q) (with 6 > 0) 
vector field satisfying 



q is L— periodic with respect to x, 
V ■ g = in TT, 
q ■ u = on dfl , 



(2.3) 

VI < z < d, / qi dx dy = 0. 

f{x, y, u) be a nonnegative function defined in 



Concerning the nonlinearity, let / 
X [0, 1], such that 

/ > 0, / is L-periodic with respect to x, and of class C^'^iVt x [0, 1]), 
V(a:,y) G H, /(x,i/,0)= /(a;,y,l)= 0, 

3pG (0,1), V(a;,y) G H, Vl-p< s < s' < l,f{x,y,s) > f{x,y,s'), 
VsG (0,1), 3{x,y)eil such that /(a;,?/,s) > 0, 

V(x,y) gH, /:(x,y,0) ^^"^'^'"^ 
with the additional assumption 



(2.4) 



lim 



> 0, 



u 



V (x, y, s) G X (0, 1), < fix, y, s) < fiix, y, 0) x s. 

We denote by ({x, y) := /^(x, y, 0), for each (x, y) G 

The set of such nonlinearities contains two particular types of functions: 



(2.5) 



The homogeneous (KPP) type: f{x,y,u) = g{u), where g is a. C^'^ function that 
satisfies: 
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^(0) = g{l) = 0, (7 > on (0, 1), g'{0) > 0, g'{l) < and < g{s) < g'{0)s in (0, 1). 

• Another type of such nonhnearities consists of functions f{x, y, u) = h{x, y).f{u), such 
that / is of the previous type, while h lies in C^'^{Q),L -periodic with respect to x, 
and positive in Q. 

Having this periodic framework, the notions of travelhng fronts and propagation were 
extended, in [2], [3], [18], [26] [28], [29], and [34] as follows: 

Definition 2.3 Let e = {e^, ■ ■ ■ , e*^) be an arbitrarily given vector in M'^. A function u = 
u{t, X, y) is called a pulsating travelling front propagating in the direction of e with an effective 
speed c ^ 0, if u is a classical solution of 

' Ut = V- {A{x, y)\/u) + y)-\/u+ /(x, y,u), t E R, (x, y) E n, 
V ■ A \Ju{t, x,y) = 0, t E M, (x, y) E dQ, 

< ykE[[LiZ,y{t,x,y)ERxn, u{t ,x,y) = u{t,x + k,y), (2.6) 

i=l ^ 

lim u(t,x,y) =0, and lim u(t,x,y) = 1, 

x-e—*—oo x-e—>+oo 

. < n < 1, 

where the above limits hold locally in t and uniformly in y and in the directions of which 
are orthogonal to e . 

2.2 Some important known results concerning the propagation 
phenomena in a periodic framework 

Under the assumptions (2.1), (2.2), (2.3), (2.4) and (2.5) set in the previous subsection, 
Berestycki and Hamel [2] proved that: having a pre-fixed unit vector e E M"^, there exists 
c*(e) > such that pulsating travelling fronts propagating in the direction e (i.e satisfying 
(2.6)) with a speed of propagation c exist if and only if c > c*(e); moreover, the pulsating 
fronts (within a speed c > c*(e)) are increasing in the time t. The value c*(e) = ^ ^ j(e) is 
called the minimal speed of propagation in the direction of e. Other nonhnearities have been 
considered in the cases of the whole space or in the general periodic framework (see [2], 
[28], [29], [32], [33], [34], [35]). 

Having the threshold value ^ ^ j(e), our paper aims to study the limits, the asymptotic 
behaviors, and the variations of some parametric quantities. These parametric quantities in- 
volve the parametric speeds of propagation of different reaction-advection-diffusion problems 
within a diffusion factor e > 0, a reaction factor i? > 0, or a periodicity parameter L. Thus, 
it is important to have a variational characterization which shows the dependance of the 
minimal speed of propagation on the coefficients A, q and / and on the geometry of the do- 
main Vt. In this context, Berestycki, Hamel, and Nadirashvili [3] gave such a formulation for 
CQ^gj(e) involving elliptic eigenvalue problems. We recall this variational characterization 
in the following theorem: 

Theorem 2.4 (Berestycki, Hamel, and Nadirashvili [3]) Let e be a fixed unit vector 
in W^. Let e = (e, 0, . . . , 0) E . Assume that n, A and f satisfy (2.1),(2.2), (2.4), and 
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(2.5). The minimal speed c*(e) = c^^^ ^(e) of pulsating fronts solving (2.6) and propagating 
in the direction of e is given by 

c*(e) =c^,A,g,/(e) = mm'^, (2.7) 

where k{X) = kQ^e,A,qxW ^■^ ^he principal eigenvalue of the operator i^Q,e,A,q,c,A which is de- 
fined by 

Ln,e,A,q,C^'^ ■= V ■ (AV^) - 2Ae ■ AVtp +q-Vip 
+ [X^eAe - \V ■ {Ae) - \q ■ e + 

acting on the set 

E = {ip E C^(f2), if) is L-periodic with respect to x and v ■ ASJt\} = \{i>Aeip) on dQ } . 

The proof of formula (2.7) is based on methods developed in [2], [7] and [9]. These are 
techniques of sub and super-solutions, regularizing and approximations in bounded domains. 

We note that in formula (2.7), the value of the minimal speed c*(e) is given in terms of the 
direction e, the domain Q, and the coefficients A,q and /„(., .,0). Moreover, it is important 
to notice that the dependence of c*(e) on the nonlinearity / is only through the derivative 
of / with respect to m at m = 0. 

Before going further on, let us mention that formula (2.7) extends some earlier results 
about front propag ation. When n = , A = Id and f = f{u) (with f{u) < f'{0)u 
in [0,1]), formula (2.7) then reduces to the well-known KPP formula c*(e) = 2^/'(0). 
That is the value of the minimal speed of propagation of planar fronts for the homogenous 
reaction-diffusion equation: Ut — Au = f{u) in M^.^ 

The above variational characterization of the minimal speed of propagation of pulsating 
fronts in general periodic excitable media will play the main role in studying the dependence 
of the minimal speed c*(e) = CQ^^ gj(e) on the coefficients of reaction, diffusion, advection 
and on the geometry of the domain. In this context, we have: 

Theorem 2.5 (Berestycki, Hamel, Nadirashvili [3]) Under the assumptions (2.1), (2.2), 
and (2.3) on f2, A, and q, let f = f{x,y,u) [respectively g = g{x,y,u)] be a nonnegative 
nonlinearity satisfying (2.4) arid (2.5). Let e be a fixed unit vector in M*^, where 1 < d < N, 

a) If fuix,y,^) < 9iix,y,0) for all {x,y) E fi, then 

Moreover if f^{x,y,0) <, ^ g^{x,y,0) m Q, then c*^ j^gj{e) < c*^ j^gg{e). 

b) The map B i— > c^^^^j(e) is increasing in B > and 

(^h,A,q,Bfi^) ^ , 

hmsup < +00. 

B~*+oo V B 



^In fact, the uniqueness, up to multiplication by a non-zero real number, of the first eigenvalue 
function of e /d /'(o) ~ k(X)'ijj together with this particular situation, yield that the princi- 
pal eigenfunction is constant and k{X) = + f (0) for all A > 0. Therefore by (2.7), we have 



c*(e) = min A+m] =2^7^. 
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c* (e) 

Furthermore, if Vt = or if uAe = on dil, then liminf ''^'^ — > 0. 
c) 



Cf2,Ag,/(^) - IK^-^) IU + 2,/ max_C(x,y) / max_ eA{x,y)e, (2.9) 
where ||(g.e)~||oo = max_{q{x,y).e)~ and s~ = max(— s,0) for each s G M. Furthermore, 

(x,y)£n 

the equality holds in (2.9) if and only if eAe and ( are constant, g.e = V . (Ae) = in Q 
and i/.Ae = on dQ {in the case when dQ ^ 0). 

d) Assume furthermore that f = f{u) and q = in then the map (3 i— c^^^o/(^) ^'^ 
increasing in 13 > 0. 

As a corollary of (2.9), we see that limsup — ■ — < C where C is a positive 



M^+oo VM 

constant. Furthermore, part d) implies that a larger diffusion speeds up the propagation in 
the absence of the advection field. 

We mention that the existence of pulsating travelling fronts in space-time periodic media 
was proved in Nolen, Xin [23, 24], Nolen, Rudd, Xin [25] and recently in Nadin [21, 22]. In 
[22], Nadin characterized the minimal speed of propagation and he studied the influence of 
the diffusion, the amplitude of the reaction term and the drift on the characterized speed. 

After reviewing some results in the study of the KPP propagation phenomena in a pe- 
riodic framework, we pass now to announce new results concerning the limiting behavior of 
the minimal speed of propagation within a small (resp. large) diffusion and reaction coeffi- 
cients (in some particular situations of the general periodic framework) and we will study the 
minimal speed as a function of the period of the coefficients in the KPP reaction-diffusion- 
advection (or reaction-diffusion) equation in the case where Q = M^. The proofs will be 
shown in details in section 8. The announced results will be applied to find the homoge- 
nization limit of the minimal speeds of propagation. We believe that this limit might help 
to find the homogenized equation in the "KPP" periodic framework (see section 9 for more 
details) . 

3 The minimal speed within small diffusion factors or 
within large period coefficients 

In this section, our problem is a reaction-diffusion equation with absence of advection terms: 
ut= (3V- {A{x, y)Vu) + f{x, y, n), t G M, (x, y) G fi. 



u ■ A Vu{t, x,y) = 0, t e M, (x, y) G dQ, 



(3.1) 



where P > 0. 

We mention that (3.1) is a reaction-diffusion problem within a diffusion matrix j3A. Let e 
be a unit direction in M."^. Under the assumptions (2.1), (2.2), (2.4) and (2.5), for each (3 > 0, 
there corresponds a minimal speed of propagation Cq^^^ qjI^) so that a pulsating front with 
a speed c and satisfying (3.1) exists if and only if c > c^,^a,o,/(^)- 
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Referring to part c) of Theorem 2.5, one gets < c^^^g ^ (e) < for any 

(3 > 0, where Mq = max_ ({x,y) and M = max_ eA{x,y)e. 

Consequently, there exists C > and independent of /5 such that 

V/? > 0, < ^^'^^'p.f^'^^ < C. (3.2) 

c* (e) 

The inequahty (3.2) leads us to investigate the limits of ^'^-^^^ — ^ _^ g ^^d as 

/? +00. The following theorem gives the precise limit when the diffusion factor tends 
to zero. However, it will not be announced in the most general periodic setting. We will 
describe the situation before the statement of the theorem: 

The domain will be in the form = M x C M^, where u CR'^ x R^-'^~'^ (rf > 0). If 
d = 0, then is a connected, open bounded subset of M^"^. While, in the case where 
1 < d < N — 1, ci;isa (Li, . . . , Lrf)-periodic open domain of M^^^ which satisfies (2.1); 
and hence, is a (/,Li, . . . , L^)— periodic subset of that satisfies (2.1) with / > and 
arbitrary. An element of = M x will be represented as 2; = (x, y) where x G M and 
y euCW^x M^-i-'^. 

The nonlinearity / = f{x, y, u), in this section, is a KPP nonlinearity defined on x [0, 1] 
that satisfies 

' / > 0, and of class C^'^{R xcJ x [0, 1]), 
/is (/, Li, . . . , Lrf)-periodic with respect to (x, yi, . . . , yd), when d> 1, 

f is /-periodic in x, when d = 0, 

y{x,y) en = RxcJ, f{x,y,0)= f{x,y,l)= 0, 

3pe (0,1), V(x,y) eU,\fl-p< s <s' <1, f{x,y,s) > f{x,y,s'), 
Vs e (0,1), 3(x, y) G such that /(x,?/, s) > 0, 

together with the assumptions 

fLix, y, 0) depends only on y; we denote by Civ) = fL{x, y, 0), V(x, y) E il. 
V(x,y) en = Rxuj, /:(x,|/,0) = Civ) > 0, (3.4) 
V(x,y,s) G n X (0,1), < fix,y,s) < ay)s. 

Notice that f^{x,y,u) is assumed to depend only on y, but f{x,y,u) may depend on x. 

Lastly, concerning the diffusion matrix, A[x,y) = A[y) = [Aij[y))i<:ij<N is a 
(with 6 > 0) symmetric matrix field whose entries are depending only on y, and satisfying 

A is (Li, . . . , Lrf)-periodic with respect to (i/i, . . . , ya), 
30 < ai < as, G cu,V^ G M^, (3.5) 

Theorem 3.1 Let e = (1, 0, . . . , 0) G and e > 0. Let n = Rxuj C R^ satisfy the form 
described in the previous page. Under the assumptions (3.3), (3.4), and (3.5), consider the 
reaction- diffusion equation 

ut{t,x,y)= eV ■ {A{y)Vu){t,x,y) + f{x,y,u), for {t,x,y) e Rx Q, 
v-A\Ju= ^ on RxRxdu. 



(3.3) 
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Assume, furthermore, that A and f satisfy one of the following two alternatives: 



or 



Then, 



3 a > 0, e u;, A{y)e = ae, 
/J(a;,Z/,0) = Civ), for all {x,y) e fi, 



f^{x,y,0) = C is constant, 

Wy G uj, A{y)e = a{y)e, where {3.i 
y (-^ a{y) is a positive, (Li, . . . , Ld)— periodic function overU. 



lim — ■ — '-r^ = 2 . /max C w max eAe. (3.9) 

e^0+ ^ V ^ V ^ 

Before going further on, we mention that the family of domains for which Theorem 3.1 
holds is wide. An infinite cylinder M x B^N^i{yQ, R) (where R > 0, and B^N-i{yQ, R) is the 
Euclidian ball of center yo and radius R) is an archetype of such domains. In these cylinders, 
u = BfiN-i^yo, R), I is any positive real number, and d = 0. The whole space is another 
archetype of the domain Q where d = N — 1, uj = and {/, Li, . . . , L^} is any family 

of positive real numbers. 

Remark 3.2 In Theorem 3.1, the domain = M x o; is invariant in the direction of e = 
(1,0..., 0) which is parallel to Ae ( in both cases (3.7) and (3.8)). Also, the assumption 
that the entries of A do not depend on x, yields that V .{Ae) = over Q. On the other 
hand, it is easy to find a diffusion matrix A and a nonlinearity / which satisfy, together, the 
assumptions of Theorem 3.1 while one of eAe{y) and ({y) is not constant. Referring to 
part c) of Theorem 2.5, one obtains: 



e 



V£ > 0, < '^-^ < 2 /max C(?/) Jmax eAe{y). 



However, Theorem 3.1 implies that 



lim — ■ — '-^ — = 2^ /max C(y)^ /max e74e(y). 
£—0+ Je V V f 



On the other hand, if f2 = M x u; as in Theorem 3.1, A = Id and / = f{u), Theorem 2.5 
yields that c* ^ ^ (e) = 2^/E^/f'{0), for all e > 0. □ 

In the same context, one can also find the limit when the diffusion factor goes to zero, 
but in the presence of an advection field in the form of shear flows: 

Theorem 3.3 Assume that e = (1, 0, ■ ■ ■ , 0) G M^, the domain Q = M. x u has the same 
form as in Theorem 3.1, and the coefficients f and A satisfy (3.3-3.4) o.rid (3.5) respectively. 
Assume, furthermore, that for all y G cJ, there exists a{y) positive so that A{y)e = a{y)e 
in UJ. Consider, in addition, an advective shear flow q = (g^(?/), 0, . . . , 0) (y G uJ) which is 
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(Li, ■ ■ ■ , Ld)— periodic with respect to y. Assume that e is a positive parameter and consider 
the parametric reaction- advection- diffusion problem 

ut = eV ■ {A{y)Vu) + q,{y)d^u{t,x,y) + f{x,y,u), t e R, {x,y) e Q, 
u ■ AVu{t,x,y) = 0, t e R, {x,y) e dQ, 

where q ^ over R xuJ and q has a zero average. Then, 

liin c^,eA,<?,/(e) = max {- q,iy)) = max(-g.e). (3.11) 

Moreover, in the same setting together with the additional assumptions: cAe and ( are both con 

one has 

c* (e) 

lim "'^^'^^'^ = max(-q(y).e) + 2v^v^, (3.12) 
£^0+ \/e ^ 



where ({y) = Co <ind eA{y)e = a for all y ^ u. 

The situation in this result is more general than that considered in part b) of Corollary 4.5 
in [4]. In details, the coefficients A and / can be both non-constant. Meanwhile, in the result 
of [4], the coefficients considered were assumed to satisfy the alternative (3.7). 

After having the exact value of lim — '■ — , we move now to investigate the limit of 

the minimal speed of propagation, considered as a function of the period of the coefficients 
of the reaction-diffusion equation set in the whole space M^, when the periodicity parameter 
tends to +oo. By making some change in variables, we will find a link between this problem 
and Theorem 3.1: 

Theorem 3.4 Let e = (1, 0, . . . , 0) G R^ . An element z G R^ is represented as z = {x, y) G 
M X ]R^"\ Assume that f = f{x,y,u) and A = A{y) satisfy (3.3), (3.4) and (3.5) with uo = 
R'^~^, d = N — 1, and I = Li = . . . = L^^i = 1. (That is, the domain and the coefficients 
of the equation are (1, 1, ... , 1) periodic with respect to y). Assume furthermore, that A and 

f satisfy either (3.7) or (3.8). For each L > 0, and {x,y) G M^, let A^{y) = A{y) and 

X 

f]^{x, y, u) = /(^) ^) u). Consider the reaction- diffusion problem 



ut{t,x,y) = V ■ {A^Vu){t,x,y) + f^{x,y,u), {t,x,y) G 



V-(A(|)V«)(t,x,2/) + /(|,f,^), (t,x,2/)GMxM^, ^^-^^^ 



whose coefficients are (L, . . . ,L) periodic with respect to (x, y) G M^. Then, 



lim c* A a f (e) = 2l max C(y) J max e.Ae(y). (3.14) 

The above theorem gives the limit of the minimal speed of propagation in the direction 
of e = (1, 0, ■ ■ ■ , 0) as the periodicity parameter L — > +oo. The domain is the whole space 
which is (L, ■ ■ ■ , L)— periodic whatever the positive number L. However, one can find 



L 



A, , Lq, , (^) ^ hm^ c;,_ ^^^^^^ (e) 



whenever g is a shear flow advection. Namely, in the same manner that Theorem 3.1 implies 
Theorem 3.4, one can prove that Theorem 3.3 implies 
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Theorem 3.5 Let e = (1, 0, . . . , 0) G M^. Assume that f = /(x, y, u) and A = A{y) satisfy 
(3.3), (3.4) and (3.5) with u = R^-\ d = N - 1, and I = Li = . . . = Ln^i = 1. (That 
is, the domain and the coefficients of the equation are (1, 1, . . . , 1) periodic with respect to 
y in M^~^j. Assume, furthermore, that for all y G M^"^, there exists a{y) positive so that 
A{y)e = a{y)e in M^"^. Let q = (g^ (y), 0, . . . , 0) for all y G ]R^~^ such that qi ^ over 
M^~"^, q is (1, ■ ■ ■ , 1)— periodic with respect to y and qi has a zero average. Then, 

Moreover, if eAe and C, are both constant oferM^"-"^, then 

r^T ^Riv 4 a f (e) = max (-g(t/).e) + 2v^V^, (3.16) 
where Co = Civ) o.''^^ a = eA{y)e for all y G 



3 7V-1 



In the proof of Theorem 3.3 (which implies Theorem 3.5), the assumption that the 
advection q is in the form of shear flows plays an important role in reducing the elliptic 
equation involved by the variational formula (8.13) below. Namely, since q = (q'i(y), 0, ■ ■ ■ , 0) 
and since e = (1, 0, ■ ■ ■ , 0), then the terms q{x, y) -V x^y^J and q{x, y) -e (in the general elliptic 
equation) become equal to qi{y)dx4' and qi{y) respectively. As a consequence, and due the 
uniqueness of the principal eigenfunction up to multiplication by a constant, we are able to 
choose ip independent of x, and hence, obtain a symmetric elliptic operator (without drift) 
whose principal eigenvalue was given by the variational formula (8.15) below (see section 8 
for more details). 

Remark 3.6 After the above explanations, we find that the techniques used to prove Theo- 
rem 3.3 which implies 3.5, will no longer work in the presence a general periodic advection 
field satisfying (2.3). 

Concerning the influence of advection, we mention that the limit of ''^''^ — as 

B 

B +00 (in the general periodic setting) is not yet given explicitly as a function of the 
direction e and the coefficients A, q and /. For more details one can see Theorem 4.1 in 
[4]. However, the problem of front propagation in an infinite cylinder with an underly- 
ing shear fiow was widely studied in Berestycki [1], Berestycki and Nirenberg [8]. In the 
case of strong advection, assume that f2 = M x c^;, where u; is a bounded smooth subset of 
M^""^, q = (g^(?/),0, ■ ■ ■ ,0), y E uj, and / = f{u) is a (KPP) nonlinearity. It was proved, 
in Heinze [16], that 

g j(e) 

lim =7, (3.17) 

where 

7 = sup / q, (y) dy, 

tpe D Juj 

D=i^tpeH\w), y iwr^y < /'(O), and j i;^ dy = 1^. 
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4 The minimal speed within large diffusion factors or 
within small period coefficients 

After having the hmit of o /(^)/ as e ^ 0"*", and after knowing that this hmit depends 

on max ({y) and max eAe{y), we investigate now the hmit of c^ma q A^I^^ ^'^^ 

diffusion factor M tends to +00, and we try to answer this question in a situation which is 
more general than that we considered in the previous section (in the case where the diffusion 
factor was going to 0"''). That is in the presence of an advection field and in a domain Q which 
satisfies (2.1) and which may take more forms other than those of section 3. We will find that 

in the case of large diffusion, the limit will depend on -l ({x,y)dx dy := — — / ({x,y)dx dy 



ieAei ' ' ' ^ 



and -f eAe{x, y)dx dy := -— / eAe{x, y)dx dy, where C denotes the periodicity cell of the 

Jc 1^1 Jc 

domain Vl. 

Theorem 4.1 Under the assumptions (2.1) for VL, (2.3) for the advection q., (2.4) and (2.5) 
for the nonlinearity f = f{x, y, u), let e be any unit direction of M*^. Assume that the diffusion 
matrix A = A{x,y) satisfies (2.2) together with V ■ Ae = over f2, and v ■ Ae = over 
dfl. For each M > and < 7 < 1/2, consider the following reaction- advection- diffusion 
equation 

Ut = MV ■ (A(x, y)Vu) + g(x, y)-Vu + f{x, y,u), t e R, (x, y) G fi, 
u ■ A Vu{t, x,y) = 0, t e M, (x, y) G dQ. 

Then 

c* . ,.-v Ae 



jj^oo ^'^^^'^/jj'"^ — " '^Jj'eAe{x,y)dxdyJ-fc{x,y)dxdy, 



' c y c 

where C is the periodicity cell of fl. 

Remarks 4.2 

• The setting in Theorem 4.1 is more general than that in Theorem 3.1, where: Q = 
M X cj, e = (1, 0, . . . , 0), and Ae = a{y)e. Under the assumptions of Theorem 3.1, the 
domain Q is invariant in the direction of Ae, which is that of e. Consequently, if u 
denotes the outward normal on dQ = M x du, one gets u ■ Ae = a{y) u ■ e = over dQ, 

d 

while V ■ (Ae) = -—aiy) = over Q. Moreover, in Theorem 3.1, we have only reaction 
ox 

and diffusion terms. That is g = 0. Therefore, considering the setting of Theorem 3.1, 
and taking (3A as a parametric diffusion matrix, one consequently knows the limits of 

c* .(e) 

"'^^'^ as /3 ^ 0"*~ (Theorem 3.1) and as /3 — +00 (Theorem 4.1). 

• The other observation in Theorem 4.1 is that the limit does not depend on the advection 
field q. This may play an important role in drawing counter examples to answer many 
different questions. For example, the variation of the minimal speed of propagation 
with respect to the diffusion factor and with respect to diffusion matrices which are 
symmetric positive definite. 
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Another important feature, in Theorem 4.1, is that the order of M in the denominator 

of the ratio — is equal to 1/2. It is independent of 7. Consequently, 

the case where the advection is null and there is only a reaction-diffusion equation 
follows, in particular, from the previous theorem. That is 



lim "'^'^^^^ — = 2J j- eAe{x,y)dx dy A j- C,{x,y)dxdy. 



• The previous point leads us to conclude that the presence of an advection with a factor 

TV/A /(^) 

M"^, where < 7 < 1/2, will have no more effect on the ratio — — — ■ — ^ as soon 

VM 

as the diffusion factor M gets very large. 

As far as the limit of the minimal speed of propagation within small periodic coefficients 
in the reaction-diffusion equation is concerned, the following theorem, which mainly depends 
on Theorem 4.1, treats this problem: 

Theorem 4.3 Let Q = M^. Assume that A = A{x,y), q = q{x,y) and f = f{x,y,u) are 
(1, . . . , 1)— periodic with respect to {x,y) G M^, and that they satisfy (2.2), (2.3), (2. 4), and 
(2.5) with Li = ... = L]s[ = 1. Let e be any unit direction 0/ M^, such that V ■ 

Ae = Q over M^. For each L > 0, let yl^(a;, y) = y4( — , — ), qj^{x,y) = q{—, —), and 

ij ij ij ij 

X y 

(x, u) = /(^) ^)) where [x, y) G M^. Consider the problem 



ut{t,x,y) = V ■ {A^Vu){t,x,y) + q^ ■ Vu{t,x,y) + f^{x,y,u), {t,x,y) G M x 

= V ■ |) Vn)(t, X, y) + |) ■ Vn(t, x, y) + /(|, |, u), 

whose coefficients are (L, . . . ,L) periodic with respect to {x, y) G M^. Then, 



nN 



(4.1) 



^M-,A„ g,,/,^^) = eAeix,y)dxdy Mcix, 



y)dx dy, 



where, in this setting, C = [0, 1] x ■ ■ ■ x [0, 1] C 



The above result gives the limit in any space dimension. It depends on the assumption 
V ■ (Ae) = in M^. However, if one takes = 1, and denotes the diffusion coefficient by 
a = a{x), X G M, then the previous result holds under the assumptions that a satisfies (2.2) 
and da/dx = in M. In other words, it holds when a is a positive constant. Thus, it is be 
interesting to mention that, in the one-dimensional case, the above limit was given in [13] 
and [17] within a general diffusion coefficient (which may be not constant over M). In details, 
assume that / = f{x,u) = (Ci^) ~ '^)'^ is ^ 1-periodic (KPP) nonlinearity satisfying (2.4) 
with (2.5), and M 3 x 1— a(x) is a 1— periodic function which satisfies < ai < a(x) < a2, 
for all X G M, where ai and a2 are two positive constants. For each L > 0, consider the 
reaction-diffusion equation 

d / X du\ r X 1 

dtu{t,x) = — {a{-)—] {t,x) + C(t)-«(^'^) ^(^'^) for (t,x)GMxM. (4.2) 
ox \ L ox J L L J 
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It was derived in [13] and, formally, in [17] that 

,1^, c; (e) = 2 ydT/cW. (4.3) 

where < a >^ denotes the harmonic mean of the map x ^ a{x) over [0, 1]. 

5 The minimal speed within small or large reaction 
coefficients 

In this section, the parameter of the reaction-advection-diffusion problem is the coefficient 
B multiplied by the nonlinearity /. In fact, it follows from Theorem 1.6 in Berestycki, 
Hamel and Nadirashvili [3] (recalled via Theorem 2.5 in the present paper) that the map 
B ^ (<i,A,q,Bfi^) I remains, with the assumption u.Ae = on dQ, bounded by two positive 
constants as B gets very large. Therefore, it is interesting to find the limit of c^,A,g,B/(6) / 
as -B — > +00 even in some particular situations. Moreover, it is important to find the limit 
of the same quantity as B ^ 0^. We start with the case where B +oo and then we move 
to that where B ^ 0~^. 

Theorem 5.1 Let e = (1,0,..., 0) E and B > 0. Assume that n = R x u C 
M^, A, and f satisfy the same assumptions of Theorem 3.1. That is, f and A satisfy (3.3), 
(3.4), and (3.5), and one of the two alternatives (3.7)-(3.8). Consider the reaction- diffusion 
equation 

ut{t,x,y)= W ■ {A{y)Wu){t,x,y) + B f{x,y,u), for {t,x,y) e Rx n, , . 
u-AVu= on RxRxduj. ^'^ 



Cr2,A,0,S/(^) 



lim — ' ' ' — = 2 /max ^(y) /max eAe{y). (5.2) 



Then, 



We mention that one can find the coefficients A, and / and the domain Q of the problem 
(5.1) satisfying all the assumptions of Theorem 5.1, which are the same of Theorem 3.1, 
including one of the alternatives (3.7)-(3.8) while one of ( and cAe is not constant. Owing 
to Theorem 1.10 in [3], it follows that 

V5 > 0, c*^,A,o,Bfie) $ 2V5 /max ({y) /max eAe{y), 
which is equivalent to saying that 

Cf2,A,0,B/(^) 



^ 2 /max C(y) /max eAe{y). 



B V y'^'^ V 

Therefore, there are heterogeneous settings in which the result found in Theorem 5.1 
does not follow trivially. 

We move now to study the limit when the reaction factor B tends to 0+. However, the 
situation will be more general than that in Theorem 5.1 because it will consider reaction- 
advection-diffusion equations rather than considering reaction-diffusion equations only: 
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Theorem 5.2 Under the assumptions (2.1) for Q, (2.3) for the advection q, (2.4) and (2.5) 
for the nonlinearity f = f{x, y, u), let e be any unit direction of M"'. Assume that the diffusion 
matrix A = A{x, y) satisfies (2.2) together with V ■ = over fi, and u ■ Ae = over dQ. 
For each B > and 7 > 1/2, consider the following reaction- advection- diffusion equation 

Ut = V- iA{x, y)Vu) + B^'qix, y) ■ Vu + B /(x, y,u), t e R, (x, y) G Q, 
u ■ A\/u{t,x,y) = 0, t e R, {x,y) e dn. 

Then 



liiii^ ^ ^J:^^ — = 2^ j- eAe{x, y)dx dy -f- ({x, y)dx dy, 



where C is the periodicity cell of Q. 



Having the above result one can mark a sample of notes: 

The order of B in the denominator of the ratio ^ j^-y ^ ]^ is independent of 

7 (it is equal to 1/2). Thus, whenever the advection is null, one gets 



lim^ ^'^'^^^ ) _ 2jj-eAe{x, y)dx dy ^ j -f C{^^ y)dxdy. 



Therefore, one concludes that the presence of an advection with a factor 5"^, where 
7 > 1/2, will have no more effect on the limit of the ratio ^ Qji.^)/^^ the 
reaction factor B gets very small. 

On the other hand, it is easy to check that the assumptions in Theorem 5.2 are more 
general than those in Theorem 5.1. Consequently, once we are in the more strict setting, 
which is that of Theorem 5.1, we are able to know both limits of c^,^,o,b/(^) / as B ^ +00 
and as i? ^ 0+. 



6 Variations of the minimal speed with respect to dif- 
fusion and reaction factors and with respect to peri- 
odicity parameters 

After having studied the limits and the asymptotic behaviors of the of the functions e t— > c^^^ gjl^) / 

c^,,^^^,^^^(e)/^/M(forverylargeMandforO<7< l/2),i?^ 5^ g, B/^")/^ 
(7 ^ 1/2) and L cj^iv q f ^^^^ where L is a periodicity parameter, we move now to 
investigate the variations of these functions with respect to the diffusion and reaction factors 
and with respect the periodicity parameter L. The present section will be devoted to discuss 
and answer these questions. 

We sketch first the form of the domain. Q C is assumed to be in the form M x 
which was taken in section 3. As a review, 1] = M x C M^, where a; C M'' x ]R^~'^~^ 
{d >0). If d = 0, the subset u; is a bounded open subset of M^"^. While, in the case where 
1 < d < N — 1, cijisa (Li, . . . , Ld)-periodic open domain of R^~'^ which satisfies (2.1); 
and hence, is a (/, Li, . . . , La)— periodic subset of that satisfies (2.1) with / > 0. An 
element of = M x cij will be represented as z = {x, y) where ?/ e a; C M*^ x M^"-'^"'^. With a 
domain of such form, we have: 
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Theorem 6.1 Let e = (1,0,..., 0) G M^. Assume that Q has the form M x which is 
described above, and that the diffusion matrix A = A{y) satisfies (3.5) together with the 
assumption 

A{x, y)e = A{y)e = a{y)e, for all (x, G M x ZJ; (6.1) 

where y ^ a{y) is a positive (Li, . . . , Ld)— periodic function defined overuJ. The nonlinearity 
f is assumed to satisfy (3.3) and (3.4). Moreover, one assumes that, at least, one of e ■ Ae 
and C is not constant. Besides, the advection field q (when it exists) is in the form q{x, y) = 
(q^iy), 0, . . . , 0) where q^ has a zero average over C, the periodicity cell of uj. For each /5 > 
consider the reaction- advection- diffusion problem 

ut = [3V ■{A{y)Vu) + ^q,{y)d^u + f{x,y,u), t G M, {x,y) G M x cu, 
V ■ A Vu{t, x,y) = 0, t e R, (x, y) G dQ. 

c* (e) 

Then the map (3 ^ ^./^^.v^gJ — decreasing in P > 0, and by Theorem 4-1, one has 



lim ""'"^'-1"'^^"^ = 2JieAe{y)dyJ-fay)dy, 



where C is the periodicity cell of uj. 

Remark 6.2 In the same setting of Theorem 6.1 but with no advection, that is = 0, we 
c* (e) 

still have (3 ^ ^'P^M — ^ decreasing map in (3 > 0. Moreover, if one of the alternatives 

VP 

(3. 7)- (3. 8) holds and there is no advection. Theorem 3.1 yields that 



lim — '■ — '-r^ = 2 /maxeAe(?/) . /maxC(y). 

/3^o+ ^J|3 V y'^'^ V y^'^ 

The preceding resuh yields another one concerned in the variation of the minimal speeds 

with respect to the periodicity parameter L. In the following, the domain will be the whole 

space M^. We choose the diffusion matrix A{x,y) = A{y), the shear flow q and reaction 

term / to be (1, . . . , l)-periodic and to satisfy some restrictions. For each L > 0, we assign 

x y X y 

the diffusion matrix Al(x, ?/) = ^(^'^)' advection field qj^{x,y) = ^(^' ^) 

non,_., - fqiu^ we a. .o.. to s.d. ™a«„. ..pect to 
periodicity parameter L, of the minimal speed cj^jv ^ /l^^)' ^^ich corresponds to the 
reaction-advection-diffusion equation within the (L, ■ ■ ■ , L)— periodic coefficients A^, q^ and 

Theorem 6.3 Let e = (1,0,..., 0) G M^. An element z e is represented as z = 
{x,y) G M X R^-^. Assume that A{x,y) = A{y) (for all {x,y) G R^ ) and f{x,y,u) sat- 
isfy ( (3.3), (3.4) and 3.5) with uj = R^-\ d = N - 1, and I = Li = . . . = L^^i = 1. 
Assume furthermore, that for all y G M^"""^, A{x,y)e = A{y)e = a{y)e, where y a{y) 
is a positive (1, . . . , l)-periodic function defined over R^~^ and that, at least, one of e ■ Ae 
and C is not constant. Let q be an advection field satisfying (2.3) and having the form 
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g(x, = . . . , 0) for each {x,y) G M^. Consider the reaction- advection- diffusion 

problem, 



(6.2) 



V(t,a;,?/) G M X M^, 
ut{t,x,y) = V ■ {A^{y)Vu){t,x,y) + {q,)^{y)d^u{t,x,y) + f^{x,y,u), 

whose coefficients are {L, . . . , L)— periodic with respect to (x, y) G M^. 
Then, the map L i— > c^n a a f i^) increasing in L > 0. 

Remark 6.4 The assumptions of Theorem 6.3 can not he fulfilled whenever N = 1. How- 
ever, assuming that N = 1 and that the function 

C ^ <a>H 



<C> 



is not identically equal to 2 (where a{x) is the diffusion factor, < a >h and < ( >a are, 
respectively, the harmonic mean of x ^ a{x) and arithmetic mean of x ^ ({x) over [0, 1]), 
it was proved, in [13] , that L c^iv „ g f (e) is increasing in L when L is close to 0. In 
particular, if a is constant and ( is not constant, or if fi is constant and a is not constant, 
then L cj^^ ^ <? / ^'^ increasing when L is close to 0. 

Concerning now the variation with respect to the reaction factor B, we have the following: 

Theorem 6.5 Assume that Q = x u and the coefficients A, q and f satisfy the same 
assumptions of Theorem 6.1. Let e = (1, . . . , 0) and for each B > 0, consider the reaction- 
advection- diffusion problem 



ut = V ■{A{y)Vu) + VBq^{y)d^u + Bf{x,y,u), t G M, (x, y) G M x 
u- AVu{t,x,y) = 0, t e R, {x,y) e dn. 

I hen, the map B t-^ is increasing m B > 0. 



B 

As a first note, we mention that Theorem 6.5 holds also in the case where there is no 
advection. On the other hand, Berestycki, Hamel and Nadirashvili [3] proved that the map 
^ ^ ^h,A,q,Bfi^) is increasing in _B > under the assumptions (2.1), (2.2), (2.3), (2.4), and 
(2.5) which are less strict than the assumptions considered in our present theorem. However, 

c* (e) 

the present theorem is concerned in the variation of the map B ' ^'^•^ — rather than 

^/B 

that oi B^ c*^,A,q,Bfie)- 

Remark 6.6 Owing to the same justifications given after Theorem 3.5, one concludes the 
importance of taking, in section 6, an advection in the form of shear flows. To study the 
variations of the minimal speeds as in Theorems 6.1, 6.3 and 6.5, hut in a more general 
framework (general advection fields, general diffusion, etc.), formula 2.1 remains an im- 
portant tool. However, we will no longer have variational formulations as (8.65) below. 
These problems remains open in the general periodic framework. 
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7 The minimal speed as a function of the positive def- 
inite diffusion matrix, Counter examples 

We studied the variation of the function (3 ^ — — '^^LJd. — the case where 1^ = M x cj, q 

vl3 _ 

is a shear flow of the form y) = {qi{y), 0, . . . , 0) on f2, while A and / satisfy (3.5), (6.1), 
(3.3) and (3.4). We obtained that the map (3 t—^ — 'Vpg,/ — jg decreasing with respect to 

/3 > in both cases: gi ^ or gi = over u. 

On the other hand, Berestycki, Hamel, and Nadirashvili [3] proved ( in part 2 of Theorem 
1.10) that: having any periodic domain Q C satisfying (2.1), q = and / = f{u), then 
the map j3 i— > Cq^^^ Qj(e) is increasing in /3 > 0. 

Having the two preceding results, there arise naturally the following two questions: 

• First: Do we still have the increasing behavior of the minimal speed with respect to 
the diffusion factor (3 in the presence of an advection, even if the nonlinearity is 
homogenous? 



Second: Owing to Theorem 2.4 (Theorem 1.1 in [3]), the map D : A 
where A describes the ordered family of positive definite matrices satisfying (2.2) (we 
say that A = A{x,y) < B = B{x,y) if and only if for each {x,y) G Q and for each 
we have zA{x,y)z < zB{x,y)z. Also, we say that A < B ii and only if 



z e 



3 7V 



for each {x,y) G and for each z G M , we have zA{x^y)z < zB{x,y)z.) is well 
defined (provided that f2, q and / satisfy (2.1), (2.3), (2.4) and (2.5)). We investigate 
the variation of the minimal speed of propagation with respect to that of the matrix 
of diffusion. More precisely, ii A = A{x,y) and B = B{x,y) are two positive definite 
matrices satisfying (2.2) and ii A < B, do we still have c; 



< Cn,B,g,/(e)? 



In fact and as it was mentioned above, we have: /3 ^ CQjSAofi^) increasing in 
/3 > 0. In other words, the map D restricted to the sub-family PD^ = {(3 A, f] > 0} which 
is generated by a prefixed matrix A is increasing. So the question becomes now: Does the 
previous conclusion remain true over the sub-family PDa in the presence of an advection ? 

The answer of the two preceding questions is negative in general. First, we prove that 
the answer to the second question is negative in general for matrices A and B such that 
A < B. We then prove, in section 7.2 that, actually, the answer is negative, in general, even 
when the diffusion matrices A and B are proportional. 



7.1 A counter example devoted to answer the second question 
Notation 7.1 For each real number b, let Aj, denote the N x N matrix having the form 

/I \ 

b 



A, 



■■. 

6/ 
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Proposition 7.2 Let e = (1, 0, . . . , 0) G M^, f2 = M x C M^, where u may or may not 

be bounded, and let q = (g^ (y), 0, . . . , 0) be a shear flow with a zero average where ^ on 
aJ. Assume that the nonlinearity f depend only on y. For each e > 0, consider the reaction- 
diffusion- advection problem 

ut(t,x,y) = da^xU + bAyU + q^{y)dxu(t,x,y) + f{y,u) inRxQ, 

= V-(^Vn) +qAy)d.u + fiy,u), (7.1) 
v,,{x,y) ■ Ai,Wx,yu{t,x,y) = v^{y) ■ Vyu{t,x,y) = Ofor{t,x,y) eRxRx du, 

where v^iiv) denotes the outward unit normal on duj at the point y G du ( t^^ix, y) = (0, i^ujiv)) 
is the outward unit normal on dQ at the point [x, y) ) and Ah is the matrix introduced in 
Notation 7.1. Then, 



where C is the periodicity cell of u. 

• Moreover, if ( is constant over oj ( say ( = (q), then 

A, q, /(^) = (-^^ (^)) + 2V^- 

• In particular, if f = f{u), then 

= and 
= max (-,.(,)) + 2 V7(0). 
Proof. Consider the following change of variables: 

X 

V{t,x,y) eRxRx u, v{t, x, y) = u{t, y). 
One then has: V(t, x,y) eRxRx u, 

vt{t, X, y) = ut{t, y), dxv{t, x, y) = -^d^ u{t, y) 

\ X X 

dxxv{t,x,y) = -dxxu{t,-j=,y) and /^yv{t,x,y) = Ayu{t,-^,y). 

Owing to the invariance of VL in the x— direction, we have the boundary condition: 
V (t, x,y) & Rx dfl, i^^{x, y) ■ Vx,yv{t, x, y) = 0. Consequently, the problem (7.1) is equiv- 
alent to the problem 

' y{t,x,y) eRxRxuj 

vt{t,x,y) =bdxxV + bAyV + Vbqi{y)dxv{t,x,y) + f{y,v), 
= bAx,yV + Vbqi{y)dxV + f{y,v) in M x M x a;, 
u^{x,y) -V x,yv{t,x,y) =0 for it, x, y) eRxRx duj. 
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Let a*^ ^ ^(e) denote the minimal speed of propagation corresponding to prob- 
lem (7.2). Referring to Theorem 4.1, and choosing 7 = 1/2, one gets 



On the other hand, a* r- (e) = V^c* a ^ f(e). Together with (7.3), we 

fi, h Id, Vb q, f y, J 



obtain that lim c* a ^ /(e) = 2w-f C{y)dy. 



For the limit as 6 ^ 0+, it follows from (3.12) in Theorem 3.3 that 

a* J- (e) 

^■s A. /(«) = i-ji ■"'■J^'^-^ = (-*(^')) + 

whenever C = Co (Co is a positive constant). 

In particular, if / = ][u) is a homogenous KPP nonlinearity, then C^[y) = Co = /'(O) for 
all y E oj. □ 

Conclusion: Let e = (1,0,..., 0) and f2 = M x tu. Choose / = /(m), and q = 

(g^(y), 0, ■ ■ ■ ,0) with / q^{y)dy = 0, and so that there exists 5 > satisfying 

Jc 

2^ri^ + 5 <ma_x(-g,(y)) + 2v/r(0)-(5. 

It follows, from Propositions (7.2), that there exist > and Mq > such that: 
VO<e<eo, c* . Ae) > max (-g,(y)) + 2 V/'(0) - 5 and 



Consequently, choosing e small enough and M large enough, it follows that Am > in 
the sense of the order relation on positive definite matrices; however. 



\AM,qJ^''>^''n,A,,qJ 



e . 



Therefore the answer of the second question is negative, in general, even when the non 
linearity / is homogenous. 



7.2 A counter example devoted to answer the first question 

In this subsection, we will show an example of a reaction-advection-diffusion problem whose 
diffusion matrix varies in the sub-family of positive definite matrices P-D/d = {/^ ^d, jS > 0}, 
where Id stands for the N x N identity matrix. In this example, we will apply an advection 
field which will destruct, even if the nonlinearity / is homogenous, the increasing behavior 
of the minimal speed with respect to (3 > (part (d) of Theorem 2.5). 
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The counter example 

Let e = (1, 0, . . . , 0) G M^, ^2 = M x C M^, where u may or may not be bounded, and 
let q = 0, . . . , 0) be a shear flow with a zero average where ^ on uJ. Assume that 

the nonlinearity / = f{u) is a homogenous "KPP" nonlinearity so that 



0<2v/f(0) + 5<max {-q^{y)) - 5, 

for some S > 0. 

Step 1. Using Theorem 4.1, with 7 = 1/2, we have 



(7.4) 



lim '^.MMvgfg,/^ =2v^. 



Thus, there exists Mq := Mq{6) > such that 



Step 2. We flx Mi > max(l, Mo(5)). Then, 
< c* 



(e) < 2v/f(0) + 5 



(7.5) 



Step 3. For the flxed number Mi, we also have a/ Mi q in the form of shear flows. Theorem 
3.3 yields that 

= max(-v/Migi(2/)) = \fM'\ max(-gi(?/)). 

e^0+ il,e Id, \/ Ml q, f J/Go; j/Go; 



Consequently, there exists = ^o(<^) > such that 



\/e < So, c* , (e) > J Ml 

n,eId,^/M[q, r ' ^ ^ 



max(-gi(y)) - 5 



> 0. 



(7.6) 



Step 4. Choosing < £1 ^ min(l,£:o)5 and owing to (7.4), (7.5) and (7.6), one then gets 
\ El Id, ^/Wiq, \ Ml Id, VMi q, 

with < £1 < Ml. 

This shows that the result of part 4 in Theorem 2.5 is no longer valid in the presence of 
an advection fleld, even if one chooses the nonlinearity / as / = f{u). 

Remark 7.3 To meet with the motivation done in the beginning of section 1, we mention 
that there appears two important features in the two counter examples which were announced 
in this section. In the counter example of subsection 7.1, the two matrices A = Am and 
B = A^, with M (resp. e) chosen sufficiently large (resp. sufficiently small), satisfy the 
properties A > B and CQ^^gj(e) < c}^^B,q,fi^)'^ however, they are not proportional (that is: 
there exists no real number a such that A = aB). Meanwhile, in the counter example of 
subsection 1.2, the matrices A = Id and B = e Id (e sufficiently small) are proportional and, 
in addition, they satisfy: A > B and CQ^^ ,jj(e) < CQ^^ gj(e). 
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8 Proofs of the announced results 



In this section, we are going to demonstrate the Theorems announced in sections 3, 4, 5, 
and 6. We will proceed in 4 subsections, each devoted to proving the results announced in 
a corresponding section. 



8.1 Proofs of Theorems 3.1, 3.3 and 3.4 

Proof of Theorem 3.1. Under the assumptions of Theorem 3.1, we can apply the varia- 
tional formula (2.7) of the minimal speed. Consequently, 



-Am ^ 

A>0 A 



ckeA,ojie) = min "--'^^^^'^ , (8.1) 



where kQ^e,eA,o,cW is the first eigenvalue (for each X, e > 0) of the eigenvalue problem 



L 



n,e, eA, 0, C, A ^ = ^n,e, eA, 0, C^^) ^) ^ ^' 
u ■ AVijj =0 onR X du. 



(8.2) 



and 



^n,e, eA, 0, C, A^(^' y) = ' (^(y) V^A(x, y)) -2e XAe ■ VV^(x, y) + 

[eX^eA{y)e - XeV ■ {A{y)e) + C(z/)] ^{x,y), 

for all (x, y) X io. 

Initially, the boundary condition in (8.2) is u ■ AS/ip = Xv ■ Ae on dQ = M x du; where 
z/(x, y) is the unit outward normal at {x, y) G dQ. However, f2 = M x is invariant in 
the direction of e which is that of Ae in both alternatives (3.7) and (3.8). Consequently, 
u ■ Ae = on d^l. 

We recall that for all A > 0, and for all £ > 0, we have g ^(A) > 0. Also, the 

first eigenfunction of (8.2) is positive over Q = R xuJ, and it is unique up to multiplication 
by a non zero constant. 

In our present setting, whether in (3.7) or (3.8) and due to the assumption (3.4), one 
concludes that the coefficients in ^ C A independent of x. Moreover, in both 
alternatives (3.7) and (3.8), the direction of Ae is the same of e = (1,0,--- ,0). On the 
other hand, since = M x c^;, then for each (x, y) G d^, we have v[x,y^ = (0; i^(^(y)), where 
i^ioiy) is the outward unit normal on du at y. Consequently, the first eigenfunction of (8.2) 
is independent of x and the eigenvalue problem (8.2) is reduced to 

V, eA, 0,C,X^- = ^ V ■ (^(2/) V0(t/)) + [e X^eA{y)e + ({v)] 0(l/) 

= ''n,e,eA, 0,C^^)<^ o^erw; (8.3) 
u{x,y)-A{y)V<l){y) = {0; u^{y)) ■ A{y)V<l>{y) = onRxdu, 

where = 0(y) is positive over uJ, L— periodic (since the domain u and the coefficients of 
^ne eA ( X Iv— periodic), unique up to multiplication by a constant, and belongs to 
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In the case where > 1, let C C M^"^ denote the periodicity cell of uj. Otherwise, (i = 
and one takes C = u;. In both cases, C is bounded. Multiplying the first line of (8.3) by 0, 
and integrating by parts over C, one gets 



e 



V(f)-A{y)V(j)dy- / [eX^eA{y)e + ({y)] (P\y) dy 
c 



c 



One also notes that, in this present setting, the operator i^n,e,eA,o,c,A is self-adjoint and its 
coefficients are (Li, . . . , L^)— periodic with respect {yi, . . . ,yd). Consequently, g ^(A) 

has the following variational characterization: 



ef V^-A{y)V^dy- I [e\^eA{y)e + Civ)] ¥^\y) dy 
Jc Jc 



c 



(8.5) 

In what follows, we will assume that (3.7) is the alternative that holds. That is, eAe = a 
is constant. The proof can be imitated easily whenever we assume that (3.8) holds. 

The function y i— > ({y) is continuous and (Li, . . . , L(i)— periodic over ZJ, whose periodicity 
cell C is a bounded subset of M^^^ (whether d = or d > 1). Let yo E C C ZJ such that 
max ({y) = C(z/o) (trivially, this also holds when ( is constant). Consequently, we have 

y&w 

e[ Vif - AVif- [ ieaX' + ((y))^' 
V ^ e H\C) \ {0}, > - [eaX' + C(z/o)] • 

Jc 

This yields that 

V£ > 0, VA > 0, -k^^^^^A, 0,C^^) ^ - + C(z/o)] • (8.6) 

Consequently, 

V. > 0, VA > 0, < Ac. + (8.7) 

A A 

C iy ) 

However, the function A i— >• Xae H — attains its minimum, over M+, at X{e) = 

X 



. This minimum is equal to 2A/C(yo) ^ \^aB. From (8.7), we conclude that 



^c.e, eA, 0,C(^(^)) 

m 

Finally, (2.7) implies that c^^sAflji^) — ^^'^ ^' ^ < 2\foLe\/ C{yo), or equiv- 



alently 



V. > 0, ^^^^^^ < 2v/^V^ 
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r2 6 A f\ ) 

We pass now to prove the other sense of the inequahty for hminf — '■ — . We will 

consider formula (8.5), and then organize a suitable function which leads us to a lower 

c* (e) 
bound of liminf ^'^ — . 

We have Civo) > 0. Let 6 be such that < 5 < Civo)- Thus < C(yo) - ^ < max C{y)- 

The continuity of C, over C C uJ, yields that there exists an open and bounded set [/ C C 
such that 

VyG Z7, C(yo)-5<C(2/). (8.9) 

Designate by if), a function in T>{C) (a C°°{C) function whose support is compact), with 

supp^/; C U, and / tp"^ = 1. One will have, 
Ju 

VA > 0, Ve > 0, 

-f'neeA C(^) ^ ^ / ■ A(y)Vz^rf2/ - /" [eX^eAiy)e + C(z/)] ^'(l/)c?l/ 
' ' ' ' ^ Ju Ju 

<e[ W ■ A(|/)VV'c^|/ - [eX'a + ({yo) - 6] [ ij\y) dy 
Ju Ju 

< e I «2|VV^p - [eX'^a + C(z/o) - 5], by (3.5), 

or equivalently 

^-^-^ > Xae + -rP(^)) (8.10) 

A A 

where I3{e) = CiVo) — 5 ~ e / a2\'Vip\'^. Choosing < e < j (this is possible), we 

a2 / 
Ju 

get p{e) > 0. 

The map A t— * Aae + — (3(e) attains its minimum, over M"*^, at A(£:) = \ ^^^^ . This 

A V ea 

minimum is equal to 2y/Fa^/ [3{e). 

Now, referring to formula (8.10), one gets 

For e small enough, ^.e. ^A, 0, C ^ l^fea ^/]3{e) for all A > 0. 

A 

Together with (2.7), we conclude that 



for e small enough, — ■ — '-r^ — 

- > 2v/^ya. (8.11) 



Consequently, 



lim inf "''^^''^^ ^ > lim inf 2 a/ /3(e)\/a 



2a/ C(yo) ~ (since ip is independent of e) 
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and this holds for all < 5 < CiUo)- Therefore, one can conclude that 



e— >0+ Y £ 

c* (e) 

Finally, the inequalities (8.8) and (8.12) imply that lim "''^ exists, and it is equal 

e-»0+ 



2v/«a/C(2/o) = 2 /maxeA(i/)e /maxC(?/). 



to 



We note that the same ideas of this proof can be easily applied in the case where the 
assumption (3.8) holds. In (3.8), we have C, is constant; however, eAe is not in general. 
Meanwhile the converse is true in the case (3.7). The little difference is that, in the case 
of (3.8), we choose the subset U (of the proof done above) around the point where eAe 
attains its maximum and then we continue by the same way used above. □ 

Proof of Theorem 3.3. We have 

On,A,,Ae) = nn /"--'"-^^^'^^^\ (8.13) 

where (due to the facts that g is a shear flow, e = (1, 0, ■ ■ ■ ,0) and e is an eigenvector of the 
matrix Aijj) for all V ^ ^ ^(A) is the principal eigenvalue of the problem 

^n,e,eyi,9,c,AV'(a;,l/) = /^t7,e,eyi,g,c(A)^(x,?/) over M X t^; 
u ■ AVip =0 on M X (9tu, 

with 

^ = eV ■{A{y)V4!)-2eXa{y)d,^+ qMdxi^ 

(8.14) 

+ [e X^eA{y)e — Xq^ (y) + ({y)] ip over M x cu. 

The uniqueness of the principal eigenfunction ip up to multiplication by a constant, 
yields that one can choose independent oi x. Hence, the elliptic operator LQ^e,eA,q,c,x can 
be reduced to the symmetric operator 

^ = eV- {A{y)Vi:) + [e X'eA{y)e - Xq, (y) + ({y)] ^. 

Consequently, 

VA > 0, Ve > 0, -kn,e,eA,q,cW = 





f Vip ■ A{y)V(pdy + X j 


I'^ Qiiy)^^- j 


f [X'eeA{y)e + ({y)] ^\y) dy 
c 


111111 

veHHC)\{0} 


c 



(8.15) 

Formula (8.15) yields that 

VA > 0, Ve > 0, -kn^e,eA,q,<:W > -A max (-gi(?/)) - X^emsixeA{y)e - maxC{y), 
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or equivalently 



1, M^ maxCiy) 

VA>0, V£>0, "-^.^.^^-g-Cl^J < (-gi(y)) + XemaxeA{y)e ' 



Putting A = A(e) = -\/ ^^^^^^^ > into the last inequality yields that 

e niaXyg^e-A(j/)e 

^a,e, eA, g,c(^) 



min — — ' < max {—qi{y)) + 2-\/i /maxe ■ y4(i/)e /max({y), 

A>o A s/ew Y j/e<^ y yeuj 

and hence, 

limsupc^ ^^ _f(e) < max . (8.16) 

e^0+ fec^ 

Now, we take yo E C {C is the periodicity cell of u) such that max^g^ (— = 
~?i(z/o) > (since q is periodic with respect to y, qi ^ and qi has a zero average) and we 
take 6 > such —qi{yo) — 5 > 0. It follows, from the continuity of qi, that there exists an 
open subset U G C such that yo G U and 

\/y e U, -qi{y) > max (-gi(y)) - 5. 

Let ip he a function in T>{C) with supp?/^ C f/, and / ^/;^ = 1. Referring to (8.15), it 

Ju 

follows that 



VA>0,V£>0, ^^^^^^^^^^f^^^ > -gi(yo) -5 + Aemine-^e + i /3(e), (8.17) 

A A 

where /3(e) = min^(y) — e / a2|V'?/'p > for a small enough e > (^2 > is the constant 

y^'^ Ju 

appearing in (3.5)). 

It follows from (8.17) that 

VA>0,V£>0, ' ^'^ > -gi(yo) -5 + 2v^ /mine- AeV^. 



A V 

Together with (8.13), and since 5 > is arbitrary, one gets 

liminfc^^^ ^(e) > -gi(?/o) = max(-gi(y)). (8.18) 

Finally, (8.16) and (8.18) complete the proof of (3.11). 

Similarly, one can use the above technics to prove (3.12). However, we will do the proof 
for the sake of completeness. First, one can easily check that 

VA > 0, Ve > 0, "■'^■^t"?-''' < ma_x {-qiiy)) + Av^maxeA(y)e + r- ■ 
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Putting X = \ l > into the last inequality yields that 

niaXyg^e-A(j/)e 

min j^' ^ < max {—qi{y)) + 2 /maxe ■ A{y)e /max(^(?/). 

A>0 Aye 'yGi^ V y&i^ y j/e^-J 

Having eAe and ( as constants, one then gets 

c* (e) 

limsup ^''^'^^'f — < max(-gi(i/)) + 2v^V^ (8-19) 



On the other hand, we take Uq E C ?,o that —qiiyo) = maxjj(— Also we take any 
positive number 6 so that < 6i < —qi (yo) . It follows, from the continuity of qi with respect 
to y, that there exist three subsets U 3 yo of C such that 

-qiiy) > (-gi(yo)) - 5 > for all y e U. 

Let ijj he a function in V{C) so that i/j"^ = 1 and t/j = on C \ U. For each e > 0, we 

Jc 

^fi,e, eA, \feq, Cp^ 



have 



= mm 

where (owing to the same above justifications) 

VA>0, V£>0, -/i;Q,e,eA,v/ig,c(A) = 





r Vv? ■ A(|/)V(^o?|/ + A j 




c 


¥.6/fl(C)\{0} 







\ ■ A{y)V'^dy + X j 




f [X^ea + Co] ¥^\y)dy 
c 


J 


f ^'^{y)dy 
c 



mm 



since eAe and C are constants. Having ip G H^{C) \ {0}, it follows that 

VA > 0, V£ > 0, 'Y^^'^ > v^(-gi(yo) - 5) + Aea + -/3(e), (8.21) 

A A 

where (3{e) = C,^ — e / a2|VV^|^ > for e > small enough. Thus, 
Ju 

> 0, > v^(-gi(2/o) - 5) + 2v/FV^v^- 

Since 5 was arbitrarily chosen, one the concludes that 

liminf ^£^^4^ > -gi(yo) + 2v/^v^. 
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Together with (8.19), the proof of (3.12) is complete. □ 
Proof of Theorem 3.4. Consider the change of variables 

v{t, X, y) = u{t, Lx, Ly), {t,x,y) eRxRx 

The function u satisfies (3.13) if and only if v satisfies 

vt{t, x,y) = ■ {A{y)Vv){t, x, y) + /(x, y, v) over M x M x (8.22) 
Consequently, 

> 0' A„0,/,(^) = ^^k^A,o,/(^) (8.23) 

2 

Taking e = 1/L , and applying Theorem 3.1 to problem (8.22), one then has 

lim ^= = lim ^ ™^ = 2 / max Civ) / max eA{y)e. (8.24) 

L-^+oo [Y 6^0+ y/e y veR'^-^ ^^y yeR'^-^ ^ ' 




Finally, (8.23) together with (8.24) complete the proof of Theorem 3.4. □ 

Proof of Theorem 3.5. Under the same change of variables considered in the proof of 
Theorem 3.4 above, one gets 

c;^,^^,^^^j^(e)= Lc;,__^^^^/e)and 

Taking e = — and using (3.11), then (3.15) follows. On the other hand, (3.12) implies 
(3.16) whenever eAe and C, are constant over M^^-*^. □ 

8.2 Proofs of Theorems 4.1 and 4.3 

Proof of Theorem 4.1. The proof will be divided into three steps: 

Step 1. According to Theorem 2.4, and since v ■ Ae = on dQ, the minimal speeds 
c* , j\/r7 /(e) are given by: 

> 0' ^kMA,M^q,&^ = A 

where k^^ j^j^ M'^ q C,^^) ^^"^ ip^^'^ denote the unique eigenvalue and the positive L- 
periodic eigenfunction of the problem 

MV ■ (AV^/;^'^0 - 2MAe ■ AV?A^'*^ + M^q ■ Wip^'^^ + [X'^MeAi - XM^'q ■ I + C]^^'^ 

= K,e,MA, M^g,c(^)^^'''' 
with V ■ AV^/'^'^ = on 9^. 
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For each A > and M > 0, let A' = Av^M, and let k^^ M'^ q C^^) ~ Ai(A',M). 

Consequently, 

VM>0, °"."A.M^«,/''L„i„MA',Af)_ (3 25) 

VM A'>o A 

where /i(A , M) and ip'^ are the first eigenvalue and the unique, positive L— periodic (with 
respect to x ) eigenfunction of 



+ 



, /2_ . _ A _ 
A eAe J — q ■ e + 



(8.26) 

^A,M _ ^(A ,M)7^^'*^ inf]. 



with u ■ AVip^ '^'^ = on dil. 

Owing to the uniqueness, up to multiplication by positive constants, of the first eigen- 
function of (8.26), one may assume that: 



VA'>0, VM>0, ||^^'*^||l2(c) =1. (8.27) 

Moreover, for each M > 0, min \ — ^ is attained at A^ > 0. Thus, 

a'>o a 

VM > 0. °"-'"-^^^./'^' = . (8.28) 

a'>o A' A^ 



The above characterization of c* , n/ri f{e)lvM will be used in the next steps in 
order to prove that Ihnmf c;;^^^^^ ^(e)/v^ (resp. limsup c;;^^^^ ^(e)/v^ ) is 



greater than (resp. less than) 2^ j- eAe{x, y)dxdy ^ j- ({x, y)dxdy; and hence, complete the 
proof. 

Step 2. Fix A' > and M > 0. We divide (8.26) by ^p^' '^'^ then, using the facts W.Ae = 
in Q and u ■ Ae = on dQ, we integrate by parts over the periodicity cell C. It follows from 
(2.3) and the L— periodicity of A, ( and ijj^ that 



jY^^:ll^^ + x" ! ~eAe+ f C = MA',M)|C|, (8.29) 
Jc U^'^^) Jc Jc 



',A ,M . y^y./,A ,M 
(^A'.Af) 

where \C\ denotes the Lebesgue measure of C. Let 



mQ = + eAe = -— f eA{x,y)edxdy and m = J- ({x,y) dx dy. 

Jc 1^ I JC Jc 



One concludes that 



VA'>0, VM>0, /i(A',M) > X'^J-eAe + J-C = >^''^rno + m, 

J n J a 
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whence 



VA > 0, WM > 0, ' > A mo + -7. 

A A 

The right side of (8.30) attains its minimum over M"*" at Ag = 




^.30) 



This minimum is 



equal to l^m^. 



.^ , ^ %MA, q, f^"^' . fi{X' , M) 

Consequently, tor any M > 0, -= = mm — > l^m^m. ihis yields 

VM A'>o A 



that 



lim inf 

M— >+oo 



M 



> 2Jj-eAe{x,y)dxdyJj-({x,y)dxdy. 



(8.31) 



Step 3. Fix a' > and M > 0. Multiply (8.26) by ip^ and integrate by parts over C. 

2 



Owing to the L— periodicity of fi, A, C, and ip ' , and due to the facts that / I ip 



in fi, and that v ■ Ae = on dQ, together with (2.3), one gets 



c 



'.A ,M 



1,V ■AB = 



M I Vip^ ■ AVtp^ + A' " / eAeiip^ ) + / C ( ^ 



/,A ,M 



C 



' Jc 



q ■ eiip 



c 



'.A ,M 



',A M 



C 



',A M 



(8.32) 



/i(A',M), 



whence 



VA' >0, VM>0, 0</i(A',M) < A'^a + /?+ — ^||(g-e) 



where a = maj^e/le(x, j/) and /3 = maj^C(^) 2/)- Together with (8.30), one gets 

(x,j/)60 (x,3/)er2 

2 2 a' 

VA' > 0, VM > 0, < A'^mo + m < /i(A',M) < A'^a + /5+ — r— II (^-e)" lU 



(8.33) 



1 A 
If 7 = then — ^|| (g • e) 

then 



A 



M 



1— II (^-e) 



^-1 



A II (9 ■ e) I loo- On the other hand, if < 7 < -, 



as M ^ +00. 



Consequently, the right side of (8.33) is bounded above by a positive constant B which 
does not depend on M and 7. This yields that 

VA' > 0, < limsup/i(A',M) < +00. 

M— >+oo 

On the other hand, it follows from (2.2) and (8.32) that VA' > 0,VM > 0, 



< 










c 




1 




< 








M 






B 




< 






M' 





;.A ,M 



C 



c 



'.A ,M 



q ■ e [ip 



/.A ,Af 
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Meanwhile, lim — =0, one then gets 



VA >0, hm / |V^^'*-'r = 0, 
VA'>0, VM>0, j(i)^^^''y = l 



(8.34) 



Fix a' > 0, and let {Mn)n be a sequence converging to +oo as n — > +oo and such 
that yu(A',M„) ^ /^''(^^") as n ^ +oo. It follows, from (8.34), that Wi^^ '^'""Wm^c) ^ 1 
as n ^ +00. Thus, the sequence (V^^ '*^")n is bounded in H^{C). Therefore, there exists 
a function ip^ G H^{C) such that, up to extraction of some subsequence, the functions 
(ijj^ '*^")„ converge in LF'{C) strong, H^(C) weak and almost everywhere in C, to the function 

,oo_ Consequently, and owing to (8.34), tp^ satisfies 



c 



1 1 

2 / r I \ 2 



|y^A,oo|2 < liminf / iW^^'^^-n =0. (8.36) 

C J Af„^+oo \Jc J 

From (8.36), it follows that for all A' > 0, the function ijj'^ '°° is almost everywhere constant 
over C. On the other hand, the elliptic regularity applied on equation (8.26) for M = Mn, 

implies that VA' > 0, the function ijj^ is continuous over C. Consequently, referring to 
(8.35), one gets 

VA' > 0, ^^'-"^ = over C. (8.37) 



Consider now equation (8.26). Fix A', take M = Mn, and integrate by parts over C. It 
follows, from (2.2), (2.3) and the assumptions V.Ae = over Q with u.Ae = on dQ, that 

/ M„ V ■ {AVip^' '^'^-) = 0, / -2A' ■ AVip^''^'^" = 0, and / q- V^^' '^'^'^ = 0. Hence, 
Jc Jc Jc 



A' 



, , q-e^jj^ + A' e-Aetjj^ + / C = /^(A', M„) / (8.38) 

" Jc Jc Jc Jc 

n 

Meanwhile, the functions "ip^ converge to the constant function 'ip'^ in LF'{C) strong; 
and hence, in L^{C) strong ( C is bounded, so L'^iC) is embedded in L^{C)). Let M„ — s> +oo 
in (8.38): 

In case 7 = 1/2, one has 



A' 



q . eifj^ = a' / q-eifj^ ^ X'lp^ / g ■ e = 



^ Jc Jc Jc 

n 

as n — >■ +CX) (from (2.3)). Also, in the case < 7 < 1/2, one trivially has 

A' 



Jc 

n 
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Moreover, eAe and ( are in L°°{C). Thus, as M„ — > +00 in (8.38), we get 



Jc Jc 



One concludes that 



VA > 0, ^ = X J- eAe + -"-^ = X mo + y. (8.39) 



/Ao,(M„) 

Whence for A = Aq = w — , one gets — -1 — = 2y/rn^. 

V "^0 Aq 
On the other hand, for all M„, 



M„ A'>o A' - Aq 



Passing M„ — >• +00, one gets limsup — ' ^-^^ < — ri — = 2y/mom, and this 



holds for all sequences {M„}„ converging to +00. Thus, 



limsup "'^^^'^^g'/ — < 2a -feAe{x,y)dxdy A Sc{x,y)dxdy. (8.41) 

Having (8.31) together with (8.41), the proof of Theorem 4.1 is complete. □ 

Proof of Theorem 4.3. We will consider the change of variables similar to that made in 
the proof of Theorem 3.4: 



v(t,x,y) = u(t, Lx, Ly), (t,x,y)&M.x 



After the same calculations done there, one gets that u satisfies (4.1) if and only if v 
satisfies 

vtit, X, y) = ^V- (A(x, y)Vv){t, x,y) + ^q- Vv{t, x, y) + /(x, y, v) over M x M^. (8.42) 
Consequently, 

> 0' = (8-43) 

On the other hand, the coefficients and the domain of problem (8.42) satisfy all the 

2 

assumptions of Theorem 4.1. Taking M = 1/L and 7 = 1/2, then (8.42) can be rewritten 
as 

vt{t,x,y) = M V ■ {A{x,y)Vv){t,x,y) + q ■ Vv{t,x,y) + f{x,y,v) over M x M^. 
In this situation, the periodicity cell of the whole space is C = [0, 1] x ■ ■ ■ x [0, 1]. 
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It follows, from Theorem 4.1, that 

^Miv 1 A i« f(^) c* 1 (e) 



L^0+ / 1 M^+oo 

VF (8-44) 

= 2 ^ j- eAe{x, y)dx dy ^ j- ({x, y)dx dy. 

Having (8.43) together with (8.44), the proof of Theorem 4.3 is complete. □ 
8.3 Proofs of Theorems 5.1 and 5.2 

Proof of Theorem 5.1. The main ideas of this proof are similar to those in the demon- 
stration of Theorem 3.1. Applying the variational formula (2.7) of the minimal speed, one 

gets 

c-„,,„,,,(e) = mm^^^^^i^, (8.45) 

' ' ' A>0 A 

where kQ^e,A,o,B(W is the first eigenvalue (for each A, 5 > 0) of the eigenvalue problem: 

Ln,e,A,o,BC,\'ip{x,y) = /c^, ,_^_o_sc(^) ^(2^' ?/) over M x ^; 
1^ ■ AV^ =0 onR X duj, 

and 

Ln,e,A,o,B<:,xHx. y) = V ■ {A{y)Vij{x, y)) - 2 XAe ■ V^pix, y) + 

[X^eA{y)e - A V ■ {A{y)e) + i?C(y)] V^(x, y), 

for each (x, y) G M x cj. 

We recall that for all A > 0, and for all S > 0, we have ^ q BC,^^^ ^ ^' ^^^o, the 
first eigenfunction of (8.46) is positive over f2 = M x tJ, and it is unique up to multiplication 
by a non zero constant. 

Moreover, whether in (3.7) or (3.8) and due to (3.4), one concludes that the coefficients 
in Lj^ ^ ^ g ^ are independent of x. Hence, the first eigenfunction of (8.46) is independent 
of X and the eigenvalue problem (8.46) is reduced to 

VaaSC, ■= V ■ + [A^eA(y)e + 5C(y)]0(y) 

= kn,e,A,o,Bci^)'P over clj; (8.47) 
v{x,y)-A{y)V(t){y) = (0; i.^(y)) ■ A(y)V0(y) = on M x 9^, 

where = (j){y) is positive over cJ, L— periodic (since the domain uj and the coefficients of 
-^Qe AO BC, A //—periodic), unique up to multiplication by a constant, and belongs to 

In the case where c? > 1, let C C ^ denote the periodicity cell of uj. Otherwise, d = 
and one takes C = tu. In both cases, C is bounded. Multiplying the first line of (8.47) by 0, 
and integrating by parts over C, one gets 



V0-A(2/)V0%- / [X^eA{y)e + i?C(y)] (f'{y)dy 

kn,e,A,o,Bc(>^) = ^ f • (8-48) 

0^(1/) dy 



c 
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One also notes that, in this present setting, the operator ^ is self-adjoint and 

its coefficients are (Li, . . . , L^)— periodic with respect {yi, . . . ,yd) - Consequently, — e a o B(^'^) 
has the following variational characterization: 



c 



V^-A{y)V^dy- [ [X^eA{y)e + B ({y)] v\y) dy 
Jc 



Jc 

(8.49) 

In what follows, we will assume that (3.7) is the alternative that holds. That is, eAe = a 
is constant. The proof can be imitated easily whenever we assume that (3.8) holds. 

The function y ({y) is continuous and (Li, . . . , L^)— periodic over uJ, whose periodicity 
cell C is a bounded subset of M^"^ (whether d = or d > 1). Let yo & C C uJ such that 
max ((y) = ({yo) (trivially, this also holds when ( is constant). Consequently, we have 

f Vif- AVif- [ iaX' + Bay)W' 
H\C)\ {0}, ^ > - [aX' + B Ciyo)] . 



f'^iy) dy 

c 



This yields that 

VB>0, VA>0, -A;^^^^^o^5^(A) > - [aX^ + B C{yo)] ■ (8.50) 

Consequently, 

V5 > 0, VA > 0, ^"-"'"'^^^^^ <Xa+ (8.51) 

A A 

However, the function A t-^ Aa + (i?C(?/o)/A) attains its minimum, over M+, at X{B) = 



BCiyo) 



This minimum is equal to 2y^B(^(y(j) x y/a. 



a 

Ke AO Bd^^^)) 

From (8.51), we conclude that: — ' S^. < 2-\/ B a\/ ({yo)- 

X[B) 

Finally, (2.7) implies that 

t N • ^n,e,A,o,BC('^) I . . 

%A,o,Bf(''^ = f^^^o A 2v^«VC(yo), 

or equivalently 

Vi? > 0, < 2v/^VC(y^ (8.52) 

V B 



il yl B f\ ) 

We pass now to prove the other sense of the inequality for liminf — ' ' ' . We will 

consider formula (8.5), and then organize a suitable function which leads us to a lower 
bound of liminf — ' ' ' . 

We have C{yo) > 0. Let 6 be such that < 6 < C{yo)- Thus < C{yo) — S < max C{y). 

The continuity of over C C uJ, yields that there exists an open and bounded set t/ C C 
such that 

C(yo)-5<C(y), Vye F. (8.53) 
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Designate by ip, a function in 'D{C) (a C°°(C) function whose support is compact), with 

supp?/' C U, and / ip"^ = 1. One will have, 
Ju 

VA > 0, V5 > 0, 

< /" ■ A{y)ViPdy - [\^a + B (C(yo) - 5)] (by (8.53)) 
Ju 

< [ a2|V^|2 - [X'a + B (C(yo) - 5)] by (3.5), 



or equivalently 



^"--''•^^W >A„ + ^p(B), (8.54) 



A - A 

B '^'''^ 



where p(-B) = C(2/o) ~ <^ — 75- / VV'I^- Choosing B large enough, we get p{B) > (this is 



possible since C(yo) — 5 > and also / a;2|V'?/'|^ > 0). The map A t— Aa + — p{B) attains 

Ju A 

its minimum, over at A(£:) = \ . This minimum is equal to 2V B a ^J~p{B). 



Now, referring to formula (8.54), one gets: 

for B large enough, -^C(-^) ^ 2^/ B a \/ p{B) for all A > 0. 

A 

Together with (2.7), we conclude that 

for B large enough, — ' — > 2\/ p{B)\foL. (8.55) 

V B 

Consequently, 



lim inf ^ > liminf 2\/ p[B\\foL 



= 2 ^/^(yjjy^^ya (since ip is independent of B), 
and this holds for all < 5 < ({yo)- Therefore, one can conclude that 

liminf > 2v^v/C(y^ (8.56) 



Finally, the inequalities (8.52) and (8.56) imply that lim — ' ' '_ exists, and it is 



equal to 2y^>y ({yo) = 2 . /max eA{y)e . /max ({y) . 

y Lu y ui 

The above proof was done while assuming that the alternative (3.7) holds. The same 
ideas of this proof can be easily applied in the case where alternative (3.8) holds. In (3.8), 
we have ( is constant; however, eAe is not in general. Meanwhile the converse is true in the 
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case (3.7). The little difference is that, in the case of (3.8), we chsose the subset U (of the 
proof done above) around the point yo where eAe attains its maximum and then we continue 
by the same way used above. □ 

Proof of Theorem 5.2. According to Theorem 2.4, and since u ■ Ae = on dQ, the 
minimal speeds c* nf{^) given by: 

where ^ B'^ q B(^^^ '^^'^ denote the unique eigenvalue and the positive L-periodic 
eigenfunction of the problem 

V ■ (AV^^'^) - 2Ae ■ AVi)^'^ + B^q ■ Wip^'^ + [A^ eAe - XB^q ■ e + B (] ^^'^ 

= ^neA B^ q BC(^) '^^'^ ^' ^^^^ ^ ' = ^ ' ^^i^^'^ = on dn. 



For each A > and B > 0, let A' = X/vB, and let ^ j^j ^ B^^-^) ~ 
Consequently, 

> 0, = mm (8.57) 

VB a'>o A'fi ' ^ ^ 

where /i(A', 5) and t/^'*' are the first eigenvalue and the unique, positive L— periodic (with 
respect to x) eigenfunction of 



+ 



V ■ {AVip^ - 2X'VBe ■ AV^^ + B^q ■ Vip^ 

(8.58) 



X'^BeAe-X'B^^^ q-B + BQ 



tfj^ 'B = ^(Y, B)tlj^ in 



with u ■ AVip^ 'B = on dn. 

Owing to the uniqueness, up to multiplication by positive constants, of the first eigen- 
function of (8.58), one may assume that: 

VA' > 0, V5 > 0, \\ij^''^\\L2^c) = 1- (8.59) 

Moreover, for each B > 0, min ' ^ is attained at A^ > 0. Thus, 

a'>o a B 

^B > 0, ^^-^^^B^^Bfje) ^ ^ ^ 

,/B A'>o X'B BX'b ^ ' 

Having the above characterization, one can now imitate the steps 2 and 3 in the proof of 
Theorem 4.1 to prove that 

liminfc* r>f(e)/VB 
B^o+ ^A,B' q, BJ^ 'I 



(resp. limsupc^^ ^^7^ j{e) / y B ) is greater than (resp. less than) 



I J eAe{x,y)dx dy ^Jcix,y)dx dy; 
and hence, complete the proof of Theorem 5.2. □ 
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8.4 Proofs of Theorems 6.1, 6.3, and 6.5 



Proof of Theorem 6.1. Referring to Theorem 2.4, it follows that for each /3 > 0, we have: 

A>0 Av^ 

where /Cf7,e,/3A, v^g,c(-^) ^^^^ eigenvalue of the problem 

Ln,e,pA,^q,(;,xHx,y) = kn,e,i3A,^q,(;W4^{x,y) over R x u; 
v.AVi) = on M X du, 



(8.61) 



where 



Ln,e,pA,^q,Q,x^ = f3V-{A{y)VtP)-2(3Xa{y)d,i;+^q,{y)d,^ 



ip over M X a;. 



The boundary condition follows so from the facts that f2 = M x cj, e = (1, 0, . . . , 0) and 
that A{y)e = a{y)e over uj. These yield that u ■ Ae = over dQ and V ■ Ae = 0. Moreover, 
for each {x,y) G dfl, we have z/(x, y) = {0;iy^{y)), where h'uiy) is the outward unit normal 
on duj at y. 

On the other hand, the function ip is positive , (Li, . . . , L^)— periodic with respect to y, 
and unique up to multiplication by non-zero constants. Meanwhile, the coefficients A, q and 
( are independent of x. Thus the eigenfunction will be independent of x and our eigenvalue 
problem is reduced to 

/5V ■ ^(Z/)) + [/? X'eA{y)e - X^q, {y) + ((y)] Hy) 

= kn,e,/3A,^q,(W'^{y) for all y e uj, (8.62) 
u{x,y) ■ A{y)V^{y) = (0; z/^(?/)) ■ A(y)V7/;(?/) = on M x ^cu. 

For each A > and /? > 0, let A' = X^/J3, and let /i;n,e,/3A, v^g,c(-^) ~ /^(-^ '/^)- Since for 

each /3 > 0, min ^'^'''^ is attained at X(/3), it follows that 

A>0 A 



^ .Am. — - 

A'>o A' 



Vp>0, — ^ = mm , (8.63) 



where yU.(A',/3) is the first eigenvalue of the problem: 

Lj^ =/?V-(A(y)V^)+ [A"eA(y)e-A'g,(y) + C(z/)l^ = /i(A',/?)V^ina., 

L -I (8.64) 

■ ylV^/' = on 9u;. 

The elliptic operator in (8.64) is self-adjoint. Consequently, the first eigenvalue 
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, P) has the following characterization:^ 
VA'>0, V/3>0, -/i(A',/3) = 



mm 



'J 


f ■ A{y)Vfdy + \' j 




f 

c 


X'\A{y)e + C{y) 




j 


c 



(8.65) 



min R{\ , 8,(fi). 

v&m(c)\{o} 



For each A and /3 > 0, i— >• i?(A , /3, </:>) attains its minimum over H^{C) \ {0} at 
the eigenfunction of the problem (8.64). On the other hand, /5 hh> R[X ,l3,ip) is increasing 
as an affine function in (3. Consequently, fixing A > and taking (3 > (3 > : 



-/i(A', /5) = i?(A', /3, ^^''/^) > i?(A', ^^''0 

> min i?(A',/5',(^) = -u(A',/3'). 

</.GHi(C)\{0} 



(8.66) 



In other words, for all A' > 0, the function (3 t— > iJ.{X' , P) is decreasing. Concerning now 
the function /5 i— > ^{e)/ \fp, one takes randomly /5 > /5' > 0, hence 



^n,/^'A,v/7../^^ _ /i(A'(/3'),/?') /i(A'(/3'),/3) 

A'(/3') ^ A'(/?') 

> min ^^^''^^ = ^^./^-^.v^g./^^^ 
-A'>o A' 



which means that the function i— > ^f^qA^I V decreasing. 

Finally, when /5 +oo, one can easily check that the hypothesis of Theorem 4.1 are 
satisfied; hence, one has the limit at +oo, and that completes the proof of Theorem 6.1. □ 

Proof of Theorem 6.3. Consider the change of variables w(t,x, y) = u{t^Lx^Ly\ for 
any (t,x,y) G M x M^. One consequently has, 

^^>0' = ^4-,J^AA,j(e). (8.67) 

Taking p = 1/l\ then 

vt (t, x,y) = p\/- iA{y)Vv) {t, x,y) + ^/p q, (y) v{t, x, y) + /(x, y, v) over M x M^. 



Owing to Theorem 6.1, the function i— > c* r- (e)/\/P is decreasing in /5 > 0. 

Besides, L 1/L^ is decreasing in L > 0. Together with (8.67), one obtains that the function 



4 f (s) is increasing in L > which completes the proof of Theorem 6.3. □ 



Proof of Theorem 6.5. Referring to Theorem 2.4, it follows that for each i? > 0, we 
have: 

J = mm — 



B A>o AV5 



^To have an idea, multiply (8.64) by the positive, (Li, . . . , Ld)— periodic function -0 and integrate by parts 
over the periodicity cell C of the the domain w. 
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Owing to the same justifications explained in the proof of Theorem 6.1, g a bc('^) 
is the first eigenvalue of the problem 

V ■ {A{y)y^{y)) + [A^e ■ Ae - XVBqM + Bdy)] ^(y) = k^,,,A,^,,B(;W^ ^ u,^^ 
u{x, y) ■ A{y)V^lj{y) = (0; uUv)) ■ A{y)Vi'{y) = on M x 9^. 



For each A > and B > 0, let A' = X/Vb and /c^g ^ y/BqB(i^) ~ /^(-^ '-^)- '^^^^ ^^^^ 
eigenvalue n{X\ B) has the following characterization: 

VA'>0,Vi?>0, - /^(^''^)- 



X'B 

Wip-A{y)Vipdy / C{y)^\y)dy 

min — — h / q^ip'^ - X' eAe(f'^ - '^'^ 



^eH\C)\{0}- Jc ' Jc A' (ggg) 

II'^IIl^c) = 1 

= min R{X , 5, if). 

^eH\C)\{Q} 
\\^\\l\c) = 1 

On the other hand, B R{X' , B, (p) is decreasing in i? > 0. Consequently, fixing A' > 
and taking < B < B , 

ipeH\C)\{0}; 
\Ml'{C) = 1 

MA', 5') 



X'B' 



In other words, for all A' > 0, the function B /i(A ,B)/XB is increasing in > 0. 
Now, we take randomly < B < B . Thus, 



— mm — 



B' A'>o X'B' X'^, X B' 

fi{X'^,,B) . MA', 5) ('n,A,VB,,Bf(^) 
> —1—2 > mm — ^ ■' 



A'^, xB - A'>o X'B /B 
which means that B c*^ ^ Bf^^^ I increasing in i? > 0. □ 



9 Applications to homogenization problems 

The reaction-advection-diffusion problem set in a heterogenous periodic domain f2 satisfying 
(2.1) generates a homogenization problem: 

Let e e M*^ be a vector of unit norm. Assume that fi, A, and / are (Li, . . . ^L^) — 
periodic and that they satisfy (2.1), (2.2), (2.3), (2.4) and (2.5). 
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For each e > 0, let fi"^ = 5 f2 and consider the following re-scales: 
V(x,?/)efi", Ae{x,y) = a(^^,^^ ,q^{x,y) = q(^^,^^ , and f^{x,y) = f (^^, 



X y 

e 



The coefficients A^, q^, and together with the domain are {e Li, . . . ,e L^)— periodic, 
and they satisfy similar properties to those of A, q, f and Q. 
Consider the parametric reaction-advection-diffusion problem 



ul{t,x,y) = V ■ {A,Vu'){t,x,y) + q^ ■ Vu' + fe{x,y,u'), t e R, {x,y)en', 
ly' -A, Vu'{t,x,y) = 0, t e M, {x,y) G dil", 



where y) denotes the outward unit normal on dQ^ at the point {x, y). 

Owing to the results found by Berestycki and Hamel in section 6 of [2], and since the 
coefficients and q^ together with the domain fi^ satisfy all the necessary assumptions, 

it follows that the problem (P^) admits a minimal speed of propagation c^e j-^{e) > 
such that (Pe) has a solution in the form of a pulsating front within a speed c if and only 
if c>c^.,^^^,,^,^^(e) >0. 

In this section, we investigate the limit of the parametric minimal speeds Cqe (e) 
(whose parameter is e) of the problems {Pe)e>o as 5 — 0"''. In other words, we search the 
limit of these minimal speeds as the periodicity cell = eC becomes a very small size. 
On the other hand, we study although not the most general setting, the variation of the map 
c^SA„g„/,(e) ine>0. 

Theorem 9.1 Let e G M*^ 6e a unit vector, and let Q C be a domain which is 

L— periodic and satisfying (2.1). Assume that A = A{x,y), q = q{x,y), and f = f{x,y,u) 
are L— periodic and that they satisfy (2.2), (2.3), (2.4) and (2.5) together with the assump- 
tions V.Ae = onQ and v.Ae = on dVt. For each e > 0, consider the problem 

ul{t,x,y) = W ■ {A,W)(t,x,y) + q, ■ W + fe{x,y,u'), t e R, {x,y)en', , . 
■ A,Vu%t,x,y) = 0, t G M, {x,y) G dn\ ^ ' 

where A^, f^ and q^ are the coefficients defined in the beginning of this section. Then, the 
minimal speed c a r (e) of pulsating travelling fronts propagating in the direction 

^ ^ 1 -^e 1 qei J £ 

of e and solving (9.1) satisfies 



J^+^hf, qeJe^^^ =2J eAe{x,y)dxdy \l -f C{x,y)dxdy, (9.2) 



c 

N 



where C is the periodicity cell of Vt and e = (e, 0, ■ ■ ■ , 0) G 

Proof. As a ffist notice, we mention that the domain Q!^ is the image of Vt by the a 
dilation whose center is the origin 0(0, ... ,0) and whose scale factor is equal to e. Conse- 
quently, 

for each e > 0, {ex, ey) G Vf if and only if (x, y) G fi, and 
{ex, ey) G dVf if and only if {x, y) G dVt. 

Moreover, 

V5 > 0, \/{x,y) G dVl, v^{ex, ey) = v{x,y). 
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Consider now, for each 5 > 0, the following change of variables 

f^(t, X, y) = u'^it, ex, ey) ; {t, x,y) eR x Q. 

One gets 

V(t, x,y) eRxQ, vfit, X, y) = ul{t, ex, ey), 
V,,y ■ {A{x, y)Vv'){t, X, y) = V,^y ■ {A,Vu'){t, ex, ey) = e^V- {A,Vu')it, ex, ey), 

and 

(ex ey\ 

i^eiex, ey) ■ [A^Vu^] {t, ex, ey) = u{x, y)- A[—, — ) Vu%t, ex, ey) 

I \e eJ ^9 3) 

= - v{x, y) ■ A{x, y)Vv'^{t, x, y) on M x dVt. 

The boundary condition in (9.1) yields that Ve{^x,ey) ■ [A^'Wu^]{t,ex,ey) = 0, for all 
{t,x,y) G M X dQ (which is equivalent to say: for all {t,ex,ey) G M x dQ'^). It follows 
from (9.3) that 

y{t,x,y) eRx dn, u ■ AVv'{t,x,y) =0. 
One can now conclude that: for each e > 0, satisfies (9.1) if and only if satisfies 

vlit, X, y) = ■ {AVv'){t, x,y) +-^q-Vv' + f{x, y, v'), t e R, {x, y) G Q, 
V ■ A Vv'{t, x,y)=0, t e R, {x, y) G dn. 

Having the assumptions (2.1), (2.2), (2.3), (2.4), and (2.5) on Q, A, q, and /, one gets 
that problem (9.4) admits, for each e > 0, a. minimal speed of propagation denoted by 
c* /in2 , (e). 

Moreover, due to the change of variables between and v^, it follows that for each 
e > 0, m"^ is a pulsating travelling front propagating in the direction of e within a speed c 

and solving (9.1) if and only if is a pulsating travelling front propagating in the direction 

c 

of e within a speed - and solving (9.4). This yields that 



where M ={l/ef. 

As e — > 0"^, the variable M — s> +00. Applying Theorem 4.1, with 7=2' ^^^^ i\iBX 



lim J}zM^22^1il — = 2.14- eAe{x, y)dx dy aI 4- ({x, y)dx dy. 



Therefore, lim Cqe a r (e) = 2* + eAe{x,y)dx dy a -I- ({x,y)dx dy, and the prooi 
of Theorem 9.1 is complete. □ 

Remark 9.2 It is worth noticing that, in formula 9.2, the homogenized speed depends on the 
averages of the diffusion and reaction coefficients, hut it does not depend on the advection. 
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We move now to study the variation of the map ^ ^ ^ f with respect 

to e > 0. In other words, we want to check the monotonicity behavior of the parametric 
minimal speed of propagation, whose parameter e > 0, as the periodicity cell of the domain 
of propagation shrinks or enlarges within a ratio e. In this study, we will consider the same 
situation of Theorem 6.1 and also the same notations introduced in the beginning of section 
9: 

Theorem 9.3 Let e = (1, . . . , 0). Assume that Q has the form M x u; where uj may or may 
not he hounded (precisely descrihed in section 3) and that the diffusion matrix A = A{y) 
satisfies (3.5) together with the assumption that e is an eigenvector of A{y) for all y & uJ, 
that is 

A{x, y)e = A{y)e = a{y)e, for all (x, y) eRxZJ; (9.6) 

where y ^-^ a{y) is a positive (Li, . . . , Ld)— periodic function defined overuJ. The nonlinearity 
f is assumed to satisfy (3.3) and (3.4). Assume further more that the advection field q (when 
it exists) is in the form q{x,y) = (q'i(y),0, . . . , 0) where has a zero average over C, the 
periodicity cell of u. For e > consider the reaction- advection- diffusion prohlem 

\Jt e M, V(x,?/) e fi" = M X eoo, 

ul{t,x,y) = V ■{A,Vu'){t,x,y) + q^ ■ W + fe{x,y,u'); (9.7) 
■ Ae Wu'it, X, y) = 0, t e M, (x, y) G dVl' . 

Then, the map e i— > c^e ^ ^ ^ {e) is increasing in e > 0. 

Proof of Theorem 9.3. For each e > 0, we consider the change of variables 

i;^(t, x, y) = u'^it, ex, ey) ; (t, x, ?/) G M x 
Owing to the justifications shown in the proof of Theorem 9.1, one consequently obtains 
^^>0' =^c; (^^,^^^^^^(e) = c*^^^^^^^ /e)/v^, (9.8) 

where /5(e) = {1/e f' . 

Applying Theorem 6.1, it follows that the map 77^: /5 ^{e)/ is decreasing 

in /5 > 0. On the other hand, the map ri^: e 1— > /3(e) is also decreasing in e > 0. Therefore, 
e Cqe a f (e), which is the composition 77^ 077^, is increasing in e > and this 

completes our proof. □ 

Other homogenization results, concerning reaction-advection-diffusion problems, were 
given in the case of a combustion-type nonlinearity / = fiu) satisfying 

3 ^ G (0, 1), /(s) = for all s G [0, d\ , f{s) > for all s G {9, 1), /(I) = 0, 
3p G (0, 1 — 6*), / is non-increasing on [1 — p, 1] . 

Consider the equation 

<(t, x) = V • {A{e-^ x)Vn^) + e-^q{e-^ x) ■ Vu' + f{u') in M^, (9.10) 

where the nonlinearity / satisfies (9.9), and the drift and diffusion coefficients q and A satisfy 
the general assumptions (2.2) and (2.3), with periodicity 1 in all variables xi, . . . ,xn- Fix a 
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unit vector e of M^. From Berestycki and Hamel [2], it follows that for each e > 0, problem 
(9.10) admits a unique pulsating front (ce,M^) such that 

u^(t, x) = 0^(x ■ e + Cst, x) 

where (j)^{s,x) is (e, e)— periodic in x that satisfies 0^(— cxo, .) = and 0^(+oo, .) = 
1. The functions are actually unique up to shifts in time, and one can assume that 
max0"(O, .) = e. 

Concerning problem (9.10), Heinze [15] proved that 

as e ^ 0~^, Cs ^ cq > 0, and u'^(t, x) uo{x ■ e + cot) weakly in H^^^, 
where (cq, uq) is the unique solution of the one-dimensional homogenized equation 
r a*no - CoUq + /(no) = inM, 

\ uo(-oo) = < no < no(+oo) = 1 in M, no(0) =6 K ■ ) 

and a* is a positive constant determined in [15]. 

In Theorem 1 of Caffarelli, Lee, Mellet [10], the homogenization limit was combined with 

the singular high activation limit for the reaction (one can also see [11] in this context) while 

the diffusion matrix was taken A = M^n. More precisely, the nonlinearity had the form 
1 u 

fe{u) = -/5(— ) with j3{s) a Lipschitz fucntion satisfying 

/3{s) > in (0, 1) and f3{s) = otherwise. 
These nonlinearities approach a Dirac mass at m = 1. 

10 Open problems 

In all the results of this paper, we deal with nonlinearities of the "KPP" type. In the 
periodic framework of this paper, pulsating travelling fronts exist also with other types of 
nonlinearities (see Theorems 1.13 and 1.14 in [2]). Namely, they exist when / = f{x,y,u) is 
of the "combustion" type satisfying: 

/ is globally Lipschitz-continuous in x M, 

V (x, y) e n, V s G (-co, 0] U [1, +oo), f{s, X, y) = 0, (10.1) 
3pG (0,1), V(x,y) Gfi, Vl-p< s <s' < l,/(x,2/,s) > /(x,2/,s'), 

and 



/ is L— periodic with respect to x, 

39 G (0,1), V(x,i/) gH, VsG [0,^], /(x,y,s) = 0, (10.2) 
Vs G {6, 1), 3 (x, y) eVL such that /(x, y, s) > 0, 

or when / = f{x,y,u) is of the "ZFK" (for Zeldovich- Frank- Kamenetskii) type satisfying 
(10.1) and 
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/ is L— periodic with respect to x, 

35 > 0, the restriction of / to H x [0, 1] is of class C^'^, (10.3) 
VsG (0,1), 3{x,y)eTl such that /(x, s) > 0. 

In particular, the "KPP" nonlinearities are of the "ZFK" type. 

Recently, El Smaily [12] gave min — max and max — min formulae for the speeds of propa- 
gation of problem (2.6) taken with a "ZFK" or a "combustion" nonlinearity. These formulae, 
together with the results of this paper, can give important estimates for the parametric min- 
imal speeds of the problem (2.6) when / is a "ZFK" nonlinearity which is not of the "KPP" 
type. Indeed, if / is a "ZFK" nonlinearity, one can find a "KPP" function h = h{x, y, u) 
such that 

V(x, y,u) eVLxR, f{x, y, u) < h{x, y, u). 
Referring to formula (1.17) in El Smaily [12], one can conclude that 

VM > 0,V5 > 0, V7 G M, c^,^^,M-g,i?/(^) ^ 'n,MA, q, Bh^^)- 
Moreover, if / is a "ZFK" nonlinearity satisfying the additional assumption 

V(x,y)en, /:(x,y,o)>o, (10.4) 

then one can find a "KPP" function g = g{x, y, u) such that (7 < / in ^2 x M, and thus 

VM > 0, Vfi > 0, V7 G M, 

As a consequence, under the assumptions that < 7 < 1/2, z/-/le = on dVL, and V-Ae = 
in r2. Theorem 4.1 implies that 



lim sup '^' l — <2A -feAe{x,y)dxdyA -fg'^{x,y,0)dxdy, (10.6) 



and 



lim inf ^^f^^' = 2 J I eAe{x, y)dx dy . -f h'^x, y, 0)dx dy > 0. (10.7) 

M-'+oo VM \jJc \jJc 

If / is a "combustion" nonlinearity, then problem (2.6) admits a solution (c, u) where 
c = ^ ^ ^. (e) > is unique and u = u(t,x,y) is increasing in t and it is unique up to a 
translation in t. Taking g as a. "KPP" nonlinearity such that g > f in Q x M. and using 
Theorem 4.1, it follows that 



lim sup ^'^"^'^^Zl'^ — < 2. -feAe{x,y)dxdy . Sg'^{x,y,0)dxdy 

VM \jJc ]jJc ^^Q^ 



together with lim inf Jli^L^J^liil > g 
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Similarly, one can get several estimates concerning the case of a small diffusion factors, 
small (resp. large) reaction factors, or small (resp. large) periodicity parameters. 

The above motivation gives several upper and lower estimates for the parametric speeds 
of propagation. However, the exact limits are not known. This leads us to ask about the 
asymptotics of the minimal speeds of propagation with respect to diffusion, reaction and 
periodicity factors in the "ZFK" case and about the asymptotics of the unique parametric 
speed of propagation in the "combustion" case. These studies should help, as it was done in 
section 9, in solving some homogenization problems in the "ZFK" case. 

Besides, Theorem 9.1 gives the limit of c*Qe a r (e) as e — > O"*". However, finding 

the homogenized equation of (9.1) in the "KPP" remains an open problem. 



11 Conclusions 

As we mentioned in the beginning of this paper, our first aim was to give a complete and 
rigorous analysis of the minimal speed of propagation of pulsating travelling fronts solving 
parametric heterogeneous reaction-advection-diffusion equations in a periodic framework. In 
the paper of Berestycki, Hamel and Nadirashvili [3], several upper and lower estimates for 
the parametric minimal speed of propagation were given (see Theorems 1.6 and 1.10 in [3]). 
However, the exact asymptotic behaviors of the minimal speed with respect to diffusion and 
reaction factors and with respect to the periodicity parameter L were not given there. In 
this paper, we determined the exact asymptotes of the minimal speed in the "KPP" periodic 
framework. In sections 3, 4 and 5, we proved that (under some assumptions on A, q, f and 
fl) the asymptotes of the parametric minimal speed are either 




2 /max C/maxeAe or 2a -J- eAe{x,y)dx dy \ -I- ({x,y)dx dy. 



(see Theorems 3.1, 3.4, 4.1, 4.3, 5.1 and 5.2 above). Moreover, we found in section 3 that 
the presence of an advection field, in the general form or in the form of shear flows, changes 
the asymptotic behavior of the minimal speed within a small diffusion (see Theorem 3.3 and 
Remark 3.6). Conversely, we proved in section 4 that the presence of a general advection 

MA q f^^^ 

field M'^q (where q satisfies (2.3)) has no effect on lim — — , . whenever < 



7 < 1/2 (see Theorem 4.1). Furthermore, we studied, in a particular periodic framework, the 
variations of the maps p i— > ^^^^ g^j^,^ ^ ^ ^ a a f ^^'^ ' " 



with respect to the positive variables (3, L and B respectively. Roughly speaking, we found 
that the first and the third maps have opposite senses of variations (see Theorems 6.1 and 
6.5). On the other hand. Theorem 6.3 and Theorem 9.3 yield that the minimal speed 
increases when the medium undergoes a dilation whose scale factor is greater than 1. 

The second aim was to find the homogenized "KPP" minimal speed. We achieved this 
goal in section 9 (Theorem 9.1) under the assumptions of free divergence on A{x,y)e and 
invariance of the domain in the direction A[x,y)e. This was an application to the results 
obtained in section 4. The found homogenized speed should play an important role in finding 
the homogenized reaction-advection-diffusion equation in the "KPP" case. In a forthcoming 
paper [13], we find also the homogenized speed in the one dimensional case but in a more 



45 



general setting (in fact, the assumption of divergence free is equivalent to the assumption 
that the diffusion term x i—* a{x) is constant over M in the case = 1). 

All the mathematical results obtained in this paper can be applied to study some spread- 
ing phenomena. Referring to the results of Weinberger [30], one can conclude that the 
spreading speed is equal to the "KPP" minimal speed of propagation in the periodic frame- 
work under some assumptions on the initial data Uq := uq{x, y) = n(0, x, y) which is defined 
on a periodic domain Q of M^. In such a setting, all our results can be applied to give 
rigorous answers on the asymptotic behavior of the parametric spreading speed with respect 
to diffusion and reaction factors and with respect to the periodicity parameter. 
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